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Preface 



This document describes the design of the accelerator and target systems for the AGS 
Super Neutrino Beam Facility. Under the direction of the Associate Laboratory Director 
Tom Kirk, BNL has established a Neutrino Working Group to explore the scientific case and 
facility requirements for a very long baseline neutrino experiment. Results of a study of the 
physics merit and detector performance was published in BNL-69395 in October 2002 [I] [2], 
where it was shown that a wide-band neutrino beam generated by a 1 MW proton beam from 
the AGS, coupled with a half megaton water Cerenkov detector located deep underground in 
the former Homestake mine in South Dakota would be able to measure the complete set of 
neutrino oscillation parameters: 

• precise determination of the oscillation parameters Anr| 2 and sin 2 2# 32 

• detection of the oscillation of — > v e and measurement of sin 2 2#i3 

• measurement of Am 21 sin2# 12 in a — > u e appearance mode, independent of the value 

of e 13 

• verification of matter enhancement and the sign of Am^ 

• determination of the CP- violation parameter Sep in the neutrino sector 

This report details the performance requirements and conceptual design of the accelerator 
and the target systems for the production of a neutrino beam by a 1.0 MW proton beam 
from the AGS. The major components of this facility include a new 1.2 GeV superconducting 
linac, ramping the AGS at 2.5 Hz, and the new target station for 1.0 MW beam. It also calls 
for moderate increase, about 30%, of the AGS intensity per pulse. Special care is taken to 
account for all sources of proton beam loss plus shielding and collimation of stray beam halo 
particles to ensure equipment reliability and personal safety. A preliminary cost estimate and 
schedule for the accelerator upgrade and target system are also included. 
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1 Introduction and Accelerator Performance 

After more than 40 years of operation, the AGS is still at the heart of the Brookhaven hadron 
accelerator complex. This system of accelerators presently comprises a 200 MeV linac for the 
pre-acceleration of high intensity and polarized protons, two Tandem Van der Graaffs for the 
pre-acceleration of heavy ion beams, a versatile Booster that allows for efficient injection of 
all three types of beams into the AGS and, most recently, the two RHIC collider rings that 
produce high luminosity heavy ion and polarized proton collisions. For several years now, the 
AGS has held the world intensity record with more than 7 x 10 13 protons accelerated in a 
single pulse 

The requirements for the proton beam for the super neutrino beam are summarized in Table 
11.11 and a layout of the upgraded AGS is shown in Figure 11.11 Since the present number of 
protons per fill is already close to the required number, the upgrade focuses on increasing 
the repetition rate and reducing beam losses (to avoid excessive shielding requirements and 
to maintain activation of the machine components at workable level). It is also important to 
preserve all the present capabilities of the AGS, in particular its role as injector to RHIC. 



Table 1.1: Performance of the 


Dresent and upgrade AGS 




Present 


Upgrade 


Average Beam Power 


0.14 MW 


1.0 MW 


Beam Energy 


24 GeV 


28 GeV 


Average Beam Current 


6 fiA 


36 //A 


Cycle Time 


2 sec 


400 ms 


Number of Protons per Fill 


7.0 x 10 13 


8.9 x 10 13 


Number of Bunches per Fill 


12 


23 


Protons per Bunch 


5.8 x 10 12 


3.72 x 10 12 


Number of Injected Turns 


190 


240 


Repetition Rate 


0.5 Hz 


2.5 Hz 


Pulse Length 


0.35 ms 


0.72 ms 


Chopping Rate 


0.75 


0.75 


Linac Average/Peak Current 


26/35 mA 


21/28 mA 



The AGS Booster was built not only to allow the injection of any species of heavy ion 
into the AGS but to allow a fourfold increase of the AGS intensity. It is one-quarter the 
circumference of the AGS with the same aperture. However, the accumulation of four Booster 
loads in the AGS takes about 0.6 sec, and is therefore not well suited for high average beam 
power operation. To minimize the injection time to about 1 msec, a 1.2 GeV linac will be used 
instead. This linac is consists the existing warm linac of 200 MeV and a new superconducting 
linac of 1.0 GeV. The multi-turn H~ injection from a source of 30 mA and 720 \i sec pulse 
width is sufficient to accumulate 9 x 10 13 particle per pulse in the AGS [I]. 

The minimum ramp time of the AGS to full energy is presently 0.5 s; this must be upgraded 
to 0.2 s to reach the required repetition rate of 2.5 Hz. The required upgrade of the AGS 
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Figure 1.1: AGS proton driver layout. 



power supply, The RF system, and other rate dependent accelerator issues will be discussed 
in Chapter 01 

The design of the target/horn configuration is shown in Figure 11.21 The material selected 
for the proton target is a Carbon-Carbon composite. It is a 3-dimensional woven material 
that exhibits extremely low thermal expansion for temperatures up to 1000°C; for higher 
temperatures it responds like graphite. This property is important for greatly reducing the 
thermo-elastic stresses induced by the beam, thereby extending the life of the target. 



l_out 




carbon-carts on 
Target Rod 



Figure 1.2: Graphite target and horn configuration. 



The target consists of an 80 cm long cylindrical rod of 12 mm diameter. The target 
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intercepts a 2 mm rms proton beam of 10 14 protons/pulse. The total energy deposited as heat 
in the target is 7.3 kJ with peak temperature rise of about 280°C. Heat will be removed from 
the target through forced convection of helium gas across its outside surface. 

The extracted proton beam uses an existing beamline at the AGS, but is then directed to 
a target station atop a constructed earthen hill. The target is followed by a downward sloping 
pion decay channel. This vertical arrangement keeps the target and decay pipe well above the 
water table in this area. The 11 degrees slope aims the neutrino beam at a water Cerenkov 
neutrino detector to be located in the Homestake mine at Lead, South Dakota. A plan view 
of the AGS facility is shown in Figure 15.21 A 3-dimensional view of the neutrino beam is 
provided in Figure 11.31 




Figure 1.3: 3-dimensional view of the neutrino beamline. The beamline is shown without 
shielding on top of the beam-line magnets and the decay tunnel. 

To assist the reading of this report and enforce consistency across chapters, the design 
parameters of each subsystem are provided in Appendix A. 

As explained in Chapter |H1 the resultant total direct cost of the 1 MW AGS Super Neutrino 
Beam Facility, not including both near and main detectors, is $218.5 M. The preliminary total 
estimated cost (TEC) is $369 M in FY03 dollars, including EDIA 15%; contingency 30%; 
BNL project overhead 13%. Escalation cannot be estimated without a project start year. 

It is estimated that three years of R & D are needed to build prototypes and complete a 
detailed engineering design that will reduce cost and improve operational reliability. This will 
be followed by 4.5 years of construction and 0.5 year of commissioning to prepare the facility 
for physics research operations. 
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2 Injector Linac 

To provide 1 MW proton beam power, the AGS has to operate at 2.5 Hz with little time 
allowed for injection. The present injector consists of the 200 MeV room temperature linac 
and 1.5 GeV Booster. It takes about 0.6 seconds to inject four Booster pulses to fill the 
AGS, which is not suitable for the upgrade operation. If a 1.2 GeV linac is used instead, the 
injection time can be reduced to less than 1 ms, allowing the AGS to cycle at the desired rate 
of 2.5 Hz. A 1.2 GeV linac injection can simultaneously fulfill the requirements of keeping 
the space charge tune shift in the AGS to be less than 0.25 and the injection losses down for 
reliable operation. 

The distance between the exit of the 200 MeV linac and the AGS injection point is about 
120 m. Only a superconducting linac (SCL) with sufficiently high gradient can meet the 
requirement of acceleration to 1.2 GeV within that distance. The superconducting linac 
technology has been used in many electron accelerators, such as LEP, CEBAF, and Tesla. 
The SNS project has successfully designed a 1.0 GeV proton SCL system with an accelerating 
gradient of about 18 MeV/m. The design of the new AGS injector linac follows closely that 
developed at SNS. This will substantially reduce the design cost and increase confidence in 
the design. 

The project described corresponds to an average SCL beam current of 37.6 fiA, that yields 
the required average beam power of 1 MW at the top energy of 28 GeV, including also a 
controlled beam loss of about 5% during multi-turn injection into the AGS. The average 
beam power in exit is 45 kW, considerably less than the 1-MW level of the equivalent 1.0- 
GeV SCL for the Spallation Neutron Source (SNS) jHj- Thus the concern about component 
activation by the induced radiation from uncontrolled beam losses is greatly reduced. The 
repetition rate of 2.5 beam pulses per second gives a beam intensity of 0.89 x 10 14 protons 
accelerated per AGS cycle; this is about 30% higher than the intensity routinely obtained 
with the present injector. At the end of an injection phase that takes about 240 turns, the 
space-charge tune depression is Av = 0.2, assuming a bunching factor (the ratio of beam peak 
current to average current), of 3. Also, with the normalized beam emittance of 100 it mm- 
mrad, the actual beam emittance at 1.2 GeV is e = 50 n mm-mrad. Obviously, the effective 
vertical acceptance of the AGS at injection is to match the final beam emittance value. The 
SCL beam pulse length is 0.72 ms, and the beam duty cycle 0.18%. 

The preliminary design of the SCL consists of three parts: (i) Low-Energy (LE), (ii) 
Medium Energy (ME), and (iii) High Energy (HE). A schematic view of the new injector is 
given in Figure 12.11 The actual location of the SCL on the BNL site is shown in Figure 12.21 
The beam leaves the present room temperature linac at the energy of 200 MeV and, after a 
bend of 17.5 degrees, enters a new 120 m long tunnel, where the SCL is located, and joins 
the AGS beamline at the location of magnet C01. The design parameters of the SCL and the 
AGS are given in Table 12.11 

2.1 Room Temperature Linac Upgrade 

The Brookhaven 200 MeV H- linac typically operates ~ 5000 hours/year in support of high 
intensity proton operation of the AGS, polarized protons for RHIC, and medical isotope 
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Figure 2.1: Layout of the 1.2-GeV Superconducting linac injector for the AGS. 

production. Some nominal operating parameters of the linac are given in Table 12.21 A 
detailed description of the linac can be found in references jB] and [7j. One can see from this 
table that the present linac can meet the requirements for beam current and repetition rate, 
and the required duty factor is less than the typical operating value. As will be discussed below, 
upgraded power supplies will be required for several systems to achieve the desired beam pulse 
width. However, since this upgrade is straightforward, the linac operation remains reliable, 
and there is room following this linac for the addition of a SCL; it was most cost effective to 
continue to use the full 200 MeV warm linac. 

The front end of the linac is shown schematically in Figure 12.31 It starts with a magnetron 
H- ion source, which produces in excess of 80 mA H~ beam at 35 keV. Following this, two 
magnetic solenoids are used to transport the H~ beam ~ lm and match it into an RFQ. The 
RFQ, operating at the linac frequency of 201.25 MHz, accelerates the beam from 35 keV to 
750 keV. After the RFQ, the beam is transported ~ 6 m to the linac. This 750 keV transport 
includes 10 quadrupoles and 3 bunchers for beam matching into the linac, and a fast beam 
chopper which allows beam chopping with ~ 10 ns rise and fall times. This chopper is a 
travelling wave structure, and the beam is chopped at a frequency to match into the Booster 
RF accelerating bucket at injection energy. A 201.25 MHz drift tube linac (DTL) accelerates 
the beam from 750 keV to 200 MeV. This linac has 9 cavities, each powered by a 5 MW 
peak power RF system. There are a total of 286 drift tubes, with a focusing electromagnetic 
quadrupole in each drift tube. 

As mentioned above, except for the beam pulse width, the present operation fulfills all 
requirements. In the future, we may consider the advantages of eliminating the long transport 
line from the RFQ to linac, in order to reduce the linac output emittance, but the present 
design assumes that this line is unchanged. Limitations in beam pulse width exist not due to 
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Table 2.1: Injector and AGS parameters for 1-MW upgrade (28 GeV). 



lIlL/ItJIll. ±jJ.11cLL, rtvc. JTUWClj IV VV 


1 .0 


Kinetic Energy, GeV 


1.2 




0.8986 


Momentum, GeV/c 


1.92 


TV /T „ x * T) ' * J * x T"! 

Magnetic Rigidity, 1-m 


0.41 


T ) 1 • j_ • T ") J T T 

Repetition Rate, Hz 


r 

2.5 


T * AT H" T"» j_ / 1 

Linac Mo. or Protons / pulse 


9.38 x 10 


Linac Duty Cycle, % 


U.l 


AGS Circumference, m 


807.076 


Revol. Frequency, MHz 


0.3338 


Revolution Period, /is 


2.996 


Bending Radius, m 


85.378 


Injection Field, kG 


0.7507 


Injection Loss, % 


5.0 


Injected Protons per Turn 


3.74 x 10 11 


Norm. Emitt., ir mm-mrad 


100 


Emittance, ir mm-mrad 


48.8 


Space-Charge Au 


0.187 



Table 2.2: Nominal operating parameters for the BNL linac. 



Output Energy 


200 MeV 


Frequency 


201.25 MHz 


Repetition Rate 


< 10Hz 


Beam Pulse Width 


< 500/is 


Nominal Duty Factor 


~ 0.3% 


Output Beam Current 


< 35 111A 


Output Emittance (rms, Normalized) 


2 7r mm mrad 


Output Energy Spread (rms) 


0.2% 
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to 20O-MeV RT Linyc 
> 



Figure 2.3: Layout of the Front-End. 

mechanical limits on the ion source or linac, but rather as a result of limits in ion source power 
supplies, RF system power supplies, and pulsed transport line and tank quadrupole power 
supplies. Therefore, in order to serve as the injector to the SCL, the following improvements 
will be required: 

1. The magnetron ion source discharge, extractor, and pulsed gas power supplies must be 
replaced by components with wider pulse capability. 

2. On the 400 kW driver stage amplifier RF systems, one must increase the 4616 plate 
capacitor bank. This entails the purchase and installation of new 35 kV capacitors in 
the existing frame. The the crowbar ignitrons and sockets must be replaced. Finally, 
we need to increase the size of the capacitors on grid power supplies. 

3. In the RF modulator system, the capacitor bank has to be increased for the 4cw25000 
anode supply, and the 8618 grid and cathode deck power supplies will need to be re- 
worked. 

4. On the 5 MW RF system, we need to replace the 60 kV capacitor banks with banks 
having more capacity; again these can be installed within existing frames. The crowbar 
ignitrons and sockets have to be replaced for 100 kA tubes. 

5. On the low level RF systems, all 400 watt solid state amplifiers will need replacement. 

6. All pulsed transport line and tank quadrupoles will have to be replaced with solid state 
units having wider pulse capabilities. 

The above upgrades are all straightforward, and we are confident that the required per- 
formance can be achieved. 
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2.2 The Superconducting Linac(SCL) 

The SCL accelerates the proton beam from 200 MeV to 1.2 GeV. The configuration and the 
design procedure of the SCL are described in detail in reference jH]. A typical sequence of 
identical periods is shown in Figure 12.41 Each period consists of a cryo-module of length 
Lcryo with an insertion of length Li ns . The insertion is needed for the placement of focusing 
quadrupoles, vacuum pumps and valves, steering magnets, beam diagnostic devices, bellows, 
and flanges. It can be either at room temperature or in a cryostat. Here we assume that the 
insertions are at room temperature. The cryo-module includes M identical cavities, each of 
N identical cells, and each having a length NL ce u, where L ce u is the length of a cell. To avoid 
coupling by the leakage of the field, cavities are separated from each other by a sufficiently 
long drift space, d. An extra drift of length L w may be added internally on both sides of the 
cryo-module to provide a transition between cold and warm regions. Thus, the length of a 
cryo-module is 

Law = MNL cM + (M — l)d + 2L W . (2.1) 

There are two symmetric intervals: a minor one, between the two middle points A and B, 
as shown in Figure 12.41 that is the interval of a cavity of length NL ce u + d; and a major one, 
between the two middle points C and D, that defines the range of a period of total length 
Lcryo + L ins . Thus, the topology of a period can be represented as a drift of length g, followed 
by M cavity intervals, and a final drift of length g, where 

g = L w + (L ins - d)/2. (2.2) 

The choice of cryo-modules with identical geometry and with the same cavity/cell config- 
uration is economical and convenient for construction. There is, nonetheless, a penalty due to 
the reduced transit-time-factors when a particle crosses cavity cells with length adjusted to a 
common central value /3q that does not correspond to the particle's instantaneous velocity. To 
minimize this affect the SCL is divided into three sections, each designed around a different 
central value /3q, and with a different cavity/cell configuration. The cell length in a section is 
fixed to be 

L C eii = AA)/2, (2.3) 

where A is the RF wavelength. We adopted an operating frequency of 805 MHz for the LE- 
section of the SCL, and 1,610 MHz for the subsequent two sections, ME and HE. The choice 
of the large RF frequency in the last two sections has been dictated by the need to achieve 
as a large accelerating gradient as possible so the SCL would fit entirely within the available 
space. The major parameters of the three sections of the SCL are given in Tables 12.31 and 

The length of the SCL depends on the average accelerating gradient. The local gradient 
has a maximum value that is limited by three causes: (1) The surface field limit, in the 
frequency range of interest between 805 and 1,610MHz, is around 40 MV/m. For a realistic 
cell shape, we set a limit on the axial electric field to 15 MV/m at 805 MHz, and 30 MV/m 
at 1,610 MHz. (2) There is a limit on the peak power provided by RF couplers that we take 
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Figure 2.4: Configuration of a proton superconducting linear accelerator. 

here not to exceed 400 kW, including a contingency of 50% to avoid saturation effects. (3) 
To make the longitudinal motion stable, we can only apply an energy gain per cryo-module 
that is a relatively small fraction of the beam energy at the exit of the cryo-module. The 
conditions for stability of motion have been derived in reference jH]. 

The proposed mode of operation is to operate each section of the SCL with the same energy 
increment. This requires the same axial field from one cryo-module to the next. To achieve 
this, and to compensate for the transit time variation from one cryo-module to next, it may 
be necessary to adjust locally the RF phase, taken here to be 30°. Also the coupling power 
may have to be adjusted according to the local transit time factor. The number of cells and 
cavities may vary from section to section, but we have found convenient here to adopt the same 
distribution in all sections. There is one klystron feeding a single coupler to a single cavity. 
The total length of the SCL injector proper from end to end is about 130 meters, including 
a 4.5-m long matching section between LE and ME sections. When averaged over the real 
estate, the actual acceleration rate is about 5 MeV/m in the LE section and 10 MeV/m in 
the ME and HE sections. Efficiencies, defined as the ratio of beam power to required total 
AC power, is relatively high for a pulsed linac, ranging between 9 and 15%. 

Negative ion stripping during transport along the SCL has been found to be very negligible, 
never exceeding a rate of 2xl0~ 8 per ion. But the final 11° bend, before injection into the 
AGS, could be a concern To control the rate of beam loss by stripping to a 10~ 5 level or 
less, the bending field should not exceed 1.25 kGauss over a total integrated bending length 
of 12 m. 

A superconducting linac is most advantageous for a continuous mode of operation (CW). 
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Table 2.3: General Parameters of the SCL. 



L/inac Section 




Mri/ 


TTTJ 


Ave. Incremental Beam Power, kW 


7.52 


15.0 


15.0 


Average Beam Current, /iA 


37.6 


37.6 


37.6 


Initial Kinetic Energy, MeV 


200 


400 


800 


Final Kinetic Energy, MeV 


400 


800 


1200 


T71 H ITT 

frequency, MHz 


805 


1610 


1610 


No. of Protons / Bunch x 10 8 


o. i U 


o. id 


O. IV 


Temperature, °K 


1 

z. 1 


1 
A. 1 


1 
A. 1 


Cells / Cavity 


8 


8 


8 


Cavities / Cryo-Module 


4 


4 


4 


Cavity Separation, cm 


32.0 


16.0 


16.0 


Cold- Warm Transition, cm 


30 


30 


30 


Cavity Internal Diameter, cm 


10 


5 


5 


Length of Warm Insertion, m 


1.079 


1.379 


1.379 


Accelerating Gradient, MeV/m 


10.5 


22.9 


22.8 


Ave. (real-estate) Gradient, MeV/m 


5.29 


9.44 


10.01 


Cavities / Klystron 


1 


1 


1 


No. of RF Couplers / Cavity 


1 


1 


1 


RF Phase Angle 


30° 


30° 


30° 


Method for Transverse Focusing 


FODO 


Doublets 


Doublets 


Betatron Phase Advance / FODO cell 


90° 


90° 


90° 


Norm, rms Emittance, ir mm-mrad 


2.0 


2.0 


2.0 


Rms Bunch Area, vr °MeV (805 MHz) 


0.5 


0.5 


0.5 
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There are two problems in the case of the pulsed- mode of operation. First, the pulsed thermal 
cycle introduces Lorentz forces that deform the cavity cells out of resonance. This can be 
controlled with a thick cavity wall strengthened on the outside by mechanical supports. The 
actual design of a cavity cell is described in detail in reference JTU]- Second, there is an 
appreciable period of time to fill the cavities with RF power before the maximum gradient 
is reached jH]. During the filling time, extra power is dissipated before the beam is injected 
into the linac. The extra amount of power required is the ratio of the filling time to the beam 
pulse length. The filling times are also shown in Table 12.41 



Table 2.4: Summary of the SCL Design. 



Linac Section 


LE 


ME 


HE 


Vplnpitv H' In 

V C±VJ> K^l \jj . kj . ±11 


W . <JVJVJ _ 


71 31 


841 8 


Out 


0.7131 


0.8418 


8986 




0.615 


0.755 


0.851 


(!p11 T iPh wf n pm 


11.45 


7 03 


7 92 


Total TVo of Ppriorlt; 


6 


9 


8 


TjPDP'tb of a FVrion 1 m 

1 J\ , 1 1 till V71 CAi J. V_^-L *_4- XIX 


6.304 


4.708 


4.994 


FODO-Cell AmDl Func Br> m 


21.52 


8.855 


8.518 


Tot p\ 1 T iPn (yf n m 

_J- w L Chi J-Jdlci 111. Ill 


37 82 


42 38 


39 96 


Connler RF Power kW (*) 

VUUUICI J. LJ- J. U VV CI ^ XV V V I j 


263 


351 


395 


Rnprffv Oriiri /Ppriod TVTpV 

J j 11V. 1. ^ — ' CAiXXX / -L Vjl J. V.7 V_X ■ _L \ _L v.y v 


33.33 


44.57 


50.10 


Total No. of Klystrons 


24 


36 


32 


Klystron Power, kW (*) 


263 


351 


395 


Z Tq, ohm/m 


378.2 


570.0 


724.2 


Qo x 10 10 


0.97 


0.57 


0.64 


Transit Time Factor, T 


0.785 


0.785 


0.785 


Ave. Axial Field, E a , MV/m 


13.4 


29.1 


29.0 


Filling Time, ms 


0.337 


0.273 


0.239 


Ave. Dissipated Power, W 


2 


11 


8 


Ave. HOM-Power, W 


0.2 


0.5 


0.4 


Ave. Cryogenic Power, W 


65 


42 


38 


Ave. Beam Power, kW 


7.52 


15 


15 


Total Ave. RF Power, kW (*) 


17 


31 


30 


Ave. AC Power for rf, kW (*) 


37 


69 


67 


Ave. AC Power for Cryo., kW 


46 


30 


27 


Total Ave. AC Power, kW (*) 


83 


99 


94 


Efficiency, % (*) 


9.05 


15.21 


16.08 



(*) Including 50% RF power contingency. 



A program was written in Visual Basic ( included with the MS Excel application), to 

calculate the beam and RF dynamics during acceleration in each of the three sections of the 

SCL. The results are displayed in Figures 12.51 to 12.101 

18 April 15, 2003 



AGS Super Neutrino Beam Facility 



Transverse Focusing 



7.fc0 



s.eo 

5.6U 



Aperture to Beam Site Ratio 




lk MlllllL-z I 



D.2E 



U.LEiC 



o:<: 



u.ulilC 



Quadrupole Gradient. 




l un K ■■. uniL 



mis BuilCh Leii-ylH (RED 3 and Spread (BLUE] 



L.fiED 



L.dLlj 



L.SED 



L.S03 



L.4ED 



L. jU j 




0.133 



C'.O-.-j 



0.053 



0.073 



O.Oi-3 



0.0E-3 



0.0^3 



I auk 



■ L-.i 



Figure 2.5: Plots (1) of behavior vs. period (tank) number of LE Section. 
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Figure 2.6: Plots (2) of behavior vs. period (tank) number of LE Section. 
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Figure 2.7: Plots (1) of behavior vs. period (tank) number of ME Section. 
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Figure 2.8: Plots (2) of behavior vs. period (tank) number of ME Section. 
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2.9: Plots (1) of behavior vs. period (tank) number of HE Section. 
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Figure 2.10: Plots (2) of behavior vs. period (tank) number of HE Section. 
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2.2.1 Transverse Focusing 

The upgrade makes use of the present 200 MeV room temperature linac, with proper mod- 
ifications as described. This linac provides a negative ion beam with an emittance of 2.0 7r 
mm-mrad. To avoid uncontrolled beam losses that may cause radiation activation, we require 
that the ratio of inner cavity radius to rms beam size is at least a factor of 6 over the length 
of the SCL. This is difficult to achieve in the ME and HE section where the inner aperture is 
of only 5 cm because of the larger RF frequency. We have thus adopted in these two sections 
transverse focusing with doublets of quadrupoles, whereas a FODO singlet sequence was found 
to be adequate in the LE section. 
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3 AGS Upgrade 

As explained in the introduction, a 1.2 GeV linac will be used for the injection into the 
AGS directly to reduce the injection time to about 1.0 msec. The results of the study of 
the multi-turn direct injection are presented in Section 13.11 The approache to increase the 
AGS repetition rate from current 0.5 Hz to 2.5 Hz is discussed in Section 13.21 In parallel 
with improvement to the power supply system is the upgrade of the RF system to raise the 
accelerating voltage from 400 kV to 1.0 MV. Finally, the eddy current effects due to increased 
repetition rate are covered in section 13.41 

All the above improvements employ proven technology and a successful implementation of 
this design can be assured. 

3.1 Multi-Turn Injection into the AGS 

The front end ion source has to operate with a 1% duty cycle at the repetition rate of 2.5 Hz 
as explained in Chapter |21 The combination of the chopper and of the RFQ prebunches the 
beam with a sufficiently small bunch length so that each beam bunch fits in the accelerating 
RF buckets of the downstream DTL, which operates at 201.25 MHz. The DTL is a room- 
temperature conventional linac that accelerates to 200 MeV. The proposed new injector for the 
AGS adds a 1.2 GeV SCL with an average output beam power of about 45 kW. The injection 
energy is still low enough to control beam losses due to stripping of the negative ions that are 
used for multiturn injection into the AGS. The duty cycle is about 0.5%. Injection into the 
AGS is modeled after the SNS scheme jS]. However, the repetition rate, and consequently the 
average beam power, is much lower here. The larger circumference of the AGS also reduces 
the number of foil traversals. Beam losses at injection into AGS are estimated to be about 3% 
controlled losses and 0.3% uncontrolled losses. This is based on a comparison with the actual 
experience in the AGS Booster and the LANL PSR and the predicted losses at the SNS using 
the quantity (Np //3 2 7 3 A), which is proportional to the Laslett tune shift, as a scaling factor. 
This is summarized in Table 13.11 As can be seen, the predicted 3% beam loss is consistent 
with both the AGS Booster and the PSR experience and the SNS prediction. 

The AGS injection parameters are summarized in Table 13.21 A relatively low RF voltage 
of 450 kV at the injection is necessary to limit the beam momentum spread at the multi-turn 
injection to be about 0.48 %, and the longitudinal emittance to be about 1.2 eVs per bunch. 
Such a small emittance is important to limit beam losses during transition crossing in the 
AGS. 

A preliminary simulation of the 360 turns injection process is shown in Figure 13.11 Without 
the second harmonic RF, some dilution of the injected particles in the phase space is inevitable. 
The bunch shape is similar to the one at the PSR of Los Alamos, with a noticeable sharp 
peak, however, a possible linac beam momentum ramping could improve this. 

The beam instability consideration is focused on two aspects. These usual for 

the AGS, at high energy, the longitudinal instability around transition, and the transverse 
instability above transition. 

The fractional beam momentum spread at transition has to be less than 0.0075 because 
of the limited momentum aperture during the transition energy jump. With the transition 
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Table 3.1: Comparison of H injection parameters. 



H Injection Parameters 


AGS Booster 


SNS 


PSR 


1 MW AGS 


Beam Power, Linac Exit (kW) 


3 


1000 


80 


45 


Kinetic Energy (MeV) 


200 


1000 


800 


1200 


No.of Protons N P (10 12 ) 


15 


100 


31 


100 


Vertical Acceptance, A (7rmm mrad) 


89 


480 


140 


55 




0.57 


6.75 


4.50 


9.56 


N P /{(3 2 -f 3 A) (10 12 /7rmm mrad) 


0.296 


0.031 


0.049 


0.190 


Total Beam Losses % 


5 


0.1 


0.3 


3 


Total Lost Beam Power W 


150 


1000 


240 


1440 


Circumference m 


202 


248 


90 


807 


Lost Beam Power per Meter W/m 


0.8 


4.0 


2.7 


1.8 



Table 3.2: AGS injection parameter. 



Injected Turns 


360 


Repetition Rate (Hz) 


2.5 


Pulse Length (ms) 


1.08 


Chopping Rate % 


0.65 


Linac Average/Peak Current (mA) 


20/30 


Momentum Spread 


±0.0015 


Norm. 95% Emittance (7r/xm) 


12 


RF Voltage (kV) 


450 


Bunch Length (ns) 


85 


Longitudinal Emittance (eVs) 


1.2 


Momentum Spread 


±0.0048 


Norm. 95% Emittance (irum) 


100 



jump, the slippage factor can be controlled to be greater than 0.002. With a bunch rms length 
of 4.25 ns and the peak current of 85 A at transition, the longitudinal impedance needs to be 
less than 11 Q to avoid longitudinal microwave instabilities. 

The measured AGS broadband impedance is about 30 Q. The broadband impedance 
mainly comes from the unshielded bellows, the vacuum chamber connections and steps and 
cavities, and also with possible contribution from BPMs and ferrite kickers. With a modest 
effort, this impedance can be reduced to be less than 10 Q, which is consistent with newly 
designed proton machines. 

In fact, if only the longitudinal microwave instability were of concern, a larger broadband 

impedance could be tolerated since the longitudinal space charge impedance of about 10 Q at 

transition, which is capacitive, has the effect of cancelling the inductive broadband impedance. 

However, the transverse instability at the high energy is more serious, even with a broadband 

impedance of 10 Q. In summary, since the intensity of 1 x 10 14 is only marginally higher 
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than the present intensity of 7 x 10 13 , the beam instability during acceleration and transition 
crossing can be avoided. 

3.2 AGS Main Power Supply Upgrade 

3.2.1 Present Mode of Operation 

The present AGS Main Magnet Power Supply (MMPS) is a fully programmable 6000 A, ±9000 
V SCR power supply. A 9 MW Motor Generator (MG), made by Siemens, is a part of the 
main magnet power supply of the accelerator. The MG permits pulsing the main magnets up 
to 50 MW peak power, while the input power of the MG itself remains constant. The highest 
power into the MG ever utilized is 7 MW, that is, the maximum average power dissipated in 
the AGS magnets has never exceeded 5 MW. 

The AGS ring comprises 240 magnets connected in series. The total resistance, R, is 0.27 Q 
and the total inductance, L, is 0.75 H. There are 12 superperiods, designed A through L, of 
20 magnets each, divided in two identical sets of 10 magnets per superperiod. 

Two stations of power supplies are each capable of delivering up to 4500 V and 6000 A. 
Every station consists of two power supplies connected in parallel. One power supply is a 12 
pulse SCR rated at +/-5000 V, 6000 A unit (P type) that is typically used for fast ramping 
during acceleration and energy recovery. The other is a lower voltage 24 pulse unit (F type), 
rated at +/-1000 V, 6000 A, that is used for flattop or slow ramping operation. The two 
stations are connected in series, with the magnet coils arranged to have a total resistance R/2 
and a total inductance of L/2. The grounding of the power supply is done only in one place, in 
the middle of station 1 or 2, through a resistive network. With this grounding configuration, 
the maximum voltage to ground in the magnets does not exceed 2500 V. The magnets are 
tested at 3 kV to ground prior to each startup of the AGS MMPS after long maintenance 
periods. 

3.2.2 Super Neutrino Beam Mode of Operation 

To cycle the AGS ring to 28 GeV at 2.5 Hz and with a ramp time of 200 ms, the magnet peak 
current is 4300 A and the peak voltage is 25 kV. Figure 13.21 displays the magnet current and 
voltage of a 2.5 Hz cycle. The total average power dissipated in the AGS magnets is estimated 
to be 3.7 MW. To limit the AGS coil voltage to ground to 2.5 kV, the AGS magnets must be 
divided into three identical sections, each powered similarly to the present AGS except that 
now the magnet loads represent only 1/6 of the total resistance and inductance. Every section 
will be powered separately with its own feed to the ring magnets and an identical system of 
power supplies, as shown in Figure 13.31 Bypass SCR's will be used across the four new P 
type stations, to bypass these units during the flattop, and ensure minimum ripple. Note that 
only station 1 will be grounded as it is done presently. 

Although the average power will not be higher than now, the peak power required is 
approximately 110 MW, exceeding the 50 MW rating of the existing MG. A new MG capable 
of providing 100 MW, may operate with 12 phases to limit, or even eliminate, the need for 
phase shifting transformers so that every power supply system would generate 24 pulses. The 
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Figure 3.2: Current and voltage cycle for 2.5 Hz operation. Also shown are the AGS dipole 
field and average power. 



generator voltage will be about 15 kV line-to-line, to limit the generator current to less than 
6000 A during pulsing. The generator will be rated at a slip frequency of 2.5 Hz. 

Running the AGS at 2.5 Hz requires that the acceleration ramp period decreases from 0.5 s 
down to 0.20 s. That is, the magnet current variation dl/dt is about 3 times larger than at 
present. Eddy current losses in the vacuum chamber are proportional to the square of (dl/dt), 
that is they are 10 times larger. However, this is still significantly below the present ramp rate 
of the AGS Booster which does not require active cooling. The increased eddy currents give 
rise to increased sextupole fields during the ramp, and will add about 20 units of chromaticity. 
The present chromaticity sextupoles will be upgraded to correct this. 

3.3 AGS RF System Upgrade 

At 2.5 Hz the peak acceleration rate is three times the present value for the AGS. With 10 
accelerating stations each station will need to supply 270 kW peak power to the beam. The 
present power amplifier design, employing a 300 kW power tetrode will be suitable to drive 
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Figure 3.3: Schematic of power supply connections to the AGS magnets for 2. 5- Hz operation. 



the cavities and supply power to the beam. The number of power amplifiers will be doubled 
so that each station will be driven by two amplifiers of the present design. This follows not so 
much from the power consideration but from the necessity to supply 2.5 times the RF voltage. 

An AGS RF station comprises four acceleration gaps surrounded by 0.35 m of ferrite 
stacks. The maximum voltage capability of a gap is not limited by the sparking threshold of 
the gap but by the ability of the ferrite to supply the magnetic induction. When the AGS 
operates with 0.5 Hz ramp the gap voltage is 10 kV. At 2.5 Hz we will need up to 25 kV per 
gap (roughly equal to the voltage from the same gap design used at the Booster, 22.5 kV) 
and this taxes the properties of the ferrite. Above a certain threshold value of B r f (20 mT 
for AGS ferrite 4L2) a ferrite becomes unstable and excessively lossy. The gap voltage at this 
B r f, max is simply given by 

V = -— / cuBrfdA = ualB rfmax In -, (3.1) 
at J a 

where u is the RF radian frequency and the variables a, b, and / are the inner and outer radius 
and length of the ferrite stack, respectively. 

The only free variable is to. If we operate the RF system at the 24th harmonic of the 
revolution frequency (9 MHz) then the required voltage of 25 kV can be achieved with a safe 
value for B r f )max of 18 mT. 

The next concern is the power dissipation in the ferrite and the thermal stress that is 

created by differential heating due to rf losses in the bulk of the material. We know from 

experience that below 300 mW/cm 3 the ferrites can be adequately cooled. The power density 
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is also proportional to B^r and is given by 



P uB 2 rf 
V ~ 2// (//Q)' 



(3.2) 



where \iQ is the quality factor of the ferrite. 

The nQ product is a characteristic of the ferrite material and depends on frequency and 
B r f. We have data on ferrite 4M2 (used in the Booster and SNS) at 9 MHz and 20 mT where 
the power dissipation is 900 mW/cm 3 . The details of the acceleration cycle determine the RF 
voltage program that is needed. For the cycle shown in Figure 13. 21 a peak voltage of 1 MV (40 
gaps each with 25 kV) is needed but for only 20 ms during acceleration. This is a duty factor 
of less than 0.05 giving an average power dissipation much less than the limit. We do not yet 
have data on the present AGS ferrite, 4L2 at 9 MHz. Characterizing 4L2 in this parameter 
regime is identified as an R&D issue but at least we know that retrofitting the AGS cavities 
with 4M2 is a fallback option. 




Figure 3.4: Bunch pattern for using harmonic 24 to create 6 bunches. 
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3.4 Eddy Current Effects 

Presently the magnet cycle of the AGS accelerator has a period of ~3.5 sec with rise time 
of ~200 msec between the injection energy and top energy. The proposed magnet cycle of 

2.5 Hz for the neutrino production operation will reduce the time between the injection and 
extraction to ~90 msec. The 75 mils thick vacuum beam pipe of the circulating beam of 
AGS (see Figure 14. 4|) is made of inconel conducting material. The increased rate of change of 
the magnetic flux (A$/dt) during the 2.5 Hz operation will generate eddy currents in both, 
the current caring coils and the vacuum chamber, that may introduce the following adverse 
effects: 

a. Ohmic heating of the coils. These Ohmic losses are added to the Ohmic losses due to the 

current that is used to excite the magnet and are due to the eddy currents generated in 
the conductors because of the rate of change of the magnetic flux (A$/dt). 

b. Ohmic heating of the chamber wall. 

c. Magnetic multipoles at the region of the circulating beam. Such magnetic multipoles that 

will vary during the acceleration of the beam may bring the circulating beam in a strong 
resonance that will results in large beam losses. 

The Ohmic losses may require cooling of the vacuum chamber wall and additional cooling of 
the magnet coil and may increase dramatically the requirements of the power supply of the 
main magnet. 

The above effects a) b) and c) have been studied by modeling the AGS magnet and 
the vacuum chamber in a 2-dimentional model, using the computer code OPERA |12j for 
electromagnetic fields calculations, and the results are reported in the following sections. 

The iron of the AGS main magnet is of silicon-steel magnetic material comprised of 1/8" 
thick laminations, which are electrically insulated and held together (sandwiched) by two 1" 
thick plates of the same magnetic material, one at each end of the magnet. The effects of 
the eddy currents generated in the iron of the magnet can also be studied by 3-dimentional 
modeling of the magnet and performing 3-dimentional electromagnetic calculations. 

3.4.1 Effect of the Eddy Currents in the Current Carrying Conductors of the 
Magnet. 

An approximate function of the excitation of the magnet (I vs time) that has been used 
in the electromagnetic calculations is shown in Figure 13.51 The minima of the function 
I(t)=Id c +Iosin(ci;t) correspond to the injection of the H~ beam and the maxima to the extrac- 
tion of the 28 GeV proton beam. The time varying magnetic flux generated by the excitation 
of the main magnet generates an EMF and therefore eddy currents on the magnets coils. The 
cross section of a coil in a main magnet is shown in Figure 14.41 and consists of a water cooled 
conductor with 32 turns (16 top and 16 bottom). 

The effect of the EMF has been calculated |T3| and can be counter balanced by increasing 
the voltage of the main magnet power supply. 
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2.5 Hz Excitation 




time | sec | 

Figure 3.5: The relative excitation strength of the main AGS magnet (I vs time), used in the 
transient electromagnetic field calculations. 



The effect of the eddy currents generated on the coil have been calculated and the results 
are shown in Figure 13.61 The Ohmic losses have been calculated for two excitations of the 
main magnet 2.5 Hz (green curve) and 5.0 Hz (red curve). The fact that the Ohmic losses 
during the I st period are comparable to those of the 3 rd period signifies that eddy currents 
acquire a steady-state very early in the magnet cycle. 

The consecutive minima of the Ohmic losses curves shown in Figure 13.61 correspond to 
the consecutive minima and maxima of the excitation curve shown in Figure 13.51 

The maxima of the Ohmic losses curves shown in Figure 13.61 correspond to the points of 
the excitation curve of Figure 13.51 where the rate of change (derivatives) is a maximum. 

As expected, the calculated integrated Ohmic losses for the 5.0 Hz operation are ~4 times 
larger than those for 2.5 Hz operation. 

In Figure 13.71 plotted are the total Ohmic losses in the coil due to both the excitation 
current and the eddy currents. Like in Figure 13 .61 only the I s * and the 3 rd periods are plotted. 

The Ohmic losses due to eddy currents are ~3% and ~12% of the Ohmic losses due to the 
excitation current, for 2.5 Hz and 5.0 Hz operation respectively. 

3.4.2 Eddy Currents in the Wall of the Vacuum Chamber 

The time varying magnetic flux generated by the excitation of the main magnet generates 
eddy currents in the wall of the vacuum chamber of the circulating beam. 

The eddy currents generated on the wall of the vacuum chamber have the following adverse 
effect: 

a. Ohmic heating on the wall of the vacuum chamber. 



b. Introduce magnetic multipoles including dipole field. 
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Figure 3.6: Ohmic losses per unit length due to eddy currents in the coil of the main magnet 
for 2.5 Hz operation (green curve) and 5.0 Hz operation (red curve). Only the Ohmic losses 
of the 1 st and the 3 rd periods are plotted in order to ascetaine that the eddy currents acquired 
a steady-state. 
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Figure 3.7: Total Ohmic losses per unit length due to eddy currents in the coil of the main 
magnet for 2.5 Hz operation (green curve) and 5.0 Hz operation (red curve). Only the Ohmic 
losses of the I s * and the 3 rd periods are plotted. 
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Figure 13.81 shows the cross section of the vacuum chamber with the magnitude of the eddy 
current density shown on the walls of the vacuum chamber with different colors. The regions 
of high current density are on the left and right edges of the vacuum chamber. The regions 
of low density are at the top and bottom of the vacuum chamber. 




Figure 3.8: The cross section of the vacuum chamber with the magnitude of the eddy current 
density shown on the walls of the vacuum chamber with different colors. The regions of high 
current density are at left and right edges of the vacuum chamber. The regions of low density 
are at the top and bottom of the vacuum chamber. 



The calculated Ohmic losses per unit length due to the eddy currents generated on the 
walls of the vacuum chamber are plotted in Figure 13.91 as a function of time. The Ohmic 
losses have been calculated for two excitations of the main magnet; 2.5 Hz (green curve) and 
5.0 Hz (red curve). The eddy current Ohmic losses as a function of time follow the same 
pattern as those in the coil of the main magnet (see previous section). 

Experimental measurements of the temperature rise of the vacuum chamber of the AGS 
have been performed for a single AGS c-type magnet when the coil of the magnet is subject 
to time varying sinusoidal current as shown in Figures I3.10| and 13.111 

3.4.3 Magnetic Multipoles Generated by the Eddy Currents 

The eddy currents generated on the wall of the vacuum chamber will affect the magnetic field 
at the circulating beam region. Due to the median plane symmetry of the magnet we only 
consider the normal multipoles which are define in this report with the relation: 
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Figure 3.9: Ohmic losses per unit length due to eddy currents in the walls of the vacuum 
chamber for 2.5 Hz operation (green curve) and 5.0 Hz operation (red curve). Only the 
Ohmic losses of the I s * and the 3 rd periods are plotted in order to ascertaine that the eddy 
currents acquired a steady-state. 



B r = b d j P sin(6')+bg W r 1 cos(26')+b se:j; r 2 cos(36 l )g6 l )+b orf r 3 cos(46 l )+b dec r 4 cos(56') + 



Figures 13.121 13. 131 and 13. 141 plot ratios b quad /b dip , b sex /b dip , and b oct /b dip , of the magnetic 
multipoles as a function of time during the 3 rd period from the start of the excitation of the 
main magnet. The calculation have shown that a steady has been achieved, therefore the 
functional form of the multipoles for any subsequent period is identical to the one shown in 
these figures. 

The 5.0 Hz excitation of the main magnet generates larger variation of the magnetic 
multipoles than the 2.5 Hz excitation. The effect of the calculated multipoles on the circulating 
beam will be studied in a separate report. 

In conclusion, the scheme for a 1 MW proton driver based on the AGS with upgraded 
injection is feasible. Indeed, the AGS beam intensity is only modestly higher than during the 
present high-intensity proton operation and, therefore, beam instability is not expected to be 
a problem during acceleration. 
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Figure 3.f0: The current I in the coil of the AGS as a function of time. 
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Figure 3.11: The temperature of the vacuum chamber as a function of time while the current 
in the coil of the magnet has the functional form shown in Figure 13.101 
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Figure 3.12: The ratio h qua d/hdi P as a function of time during the 3 rd period. The multipoles 
for any subsequent period is identical to these shown in the figure. 



b_/b di(i vs time 
2.5 Hi & iUi ^period 

1 | 1 1 1 1 1 1 1 i- 




-1j5 1 ' ' ' ' ' ' ' ' ' ' ' 1 

OA 0J6 O.S 1 1.2 14 1,6 

lime [sre] 

Figure 3.13: The ratio b sex /bdi P as a function of time during the 3 rd period. The multipoles 
for any subsequent period is identical to these shown in the figure. 
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Figure 3.14: The ratio b oc t/bdi P as a function of time during the 3 rd period. The multipoles 
for any subsequent period is identical to these shown in the figure. 
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4 AGS Injection and Extraction 

A high intensity proton beam of ~8.9xl0 13 protons at an energy of 28 GeV is required to 
irradiate the production target at a repetition rate of 2.5 Hz. This high intensity proton beam 
can be produced by injecting 1.2 GeV H~ ion from a superconducting linac (see Section 12.2)1 . 
This section provides a feasibility study of injection, acceleration and extraction of such a high 
intensity proton beam at a frequency of 2.5 HZ using the AGS machine. The following items 
have been studied here: 

1. Injection of the H~ Beam: 

a. The geometry of the H~ injection: Location of the injection region in AGS ring and 

required devices in the injection region. 

b. The calculation of the beam parameters of the circulating beam in AGS at the 

H~ injection point. These parameters are required for the matching of the beam 
parameters of the H~ injected beam to those of the circulating beam. 

2. "Stripping Foil" issues generated by the partially stripped H~ injected beam: 

3. Acceleration: 

a. The effect that the AGS-main-magnet cycle generates on the vacuum chamber of the 

AGS ring during the 2.5 Hz ramping of the main magnet. 

b. Calculations of the magnetic multipoles introduced, in the main magnetic field of the 

AGS synchrotron, by the eddy currents of the vacuum chamber. 

4. Extraction: 

a. The geometry of the extraction region. Location of the extraction region in the AGS 

ring, and required devices for beam extraction in the extraction region. 

b. The calculation of the beam parameters of the circulating beam in the AGS at the 

extraction point. These beam parameters are required for the matching of the 
beam parameters of the extracted beam, to those of the transport beam line to the 
"neutrino production" target (RTBT). 

4.1 H Injection into AGS 

Figure 12.21 is a drawing showing a sections of the AGS with the main magnets and the layout 
of the proposed superconducting linac that will accelerate the H~ ions to an energy of 1.2 
GeV for injection into the AGS. A section of the AGS-Booster ring appears in the upper left 
corner of the drawing. The H~ Injection region has been chosen to be the B20 straight section 
of the AGS. 
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4.1.1 Considerations for the 1.2 GeV H Injection into the AGS 

a. The H~ ions are to be injected from the SCJinac into the AGS using the "electron strip- 

ping" method [H]. The electrons of the injected H~ ions will be stripped by a foil 
located inside the pipe of the circulating beam of AGS. The issue is to determine the 
optimum location of the stripping foil. 

b. The integrated magnetic field {/Bdl} of the magnets comprising the High Energy Beam 

Transport (HEBT) line that transports the H~ ions from the exit of the SCL to the 
injection points (stripping-foil location) should be minimized in order to keep the ion- 
ization probability of the H6— ions jHj below lxlO -6 . For 8.9xl0 13 H~ ions injected into 
AGS at 2.5 Hz, an ionization probability of lxlO -6 is equivalent to 0.05 W of energy 
deposited into the HEBT line. This constraint on the allowed upper limit (< lxlO" 6 ) 
of the ionization probability of the H~ ions, sets an upper limit on the maximum value 
of the magnetic fields of the magnets needed to bend and focus the H~ injected beam. 

c. The HEBT line must provide space for Dipoles and Quads to be used to match the beam 

parameters of the injected H~ beam to those of the circulating proton beam of the AGS 
at the injection point. (The injection point is defined as the point where the central 
trajectory of the injected H~ beam intersects the stripping foil). 

d. In order to avoid the coherent betatron-oscillations that protons of the injected beam could 

perform while circulating in the AGS, the injected H~ beam and the central orbit of 
the circulating beam must be collinear at the injection point. This can be accomplished 
either by, locally displacing the circulating beam "locally bumping the beam" at the 
location of the injection point; and/or by adjusting the injection angle of the H~ injected 
beam. Both methods will be discussed. 

e. In order fill the phase space of the circulating beam with injected beam, the central orbit 

of the circulating beam is displaced "bumped locally" in the horizontal and vertical 
direction during the H~ beam injection. To implement the displacement of the central 
orbit of the circulating beam we use the same or similar hardware that is used to im- 
plementing the "local bump" in (e). The method to "locally bump" the beam will be 
discussed in a technical note. At time varying bump will also help reduce the number of 
"foil crossing" of the circulating beam during the time the AGS is being injected with 
H~ beam. 

f. The electrons generated by stripping of the H~ beam at the stripping foil, as well as the H° 

and H~ particles that may result from incomplete stripping, should be removed from the 
injected beam with a controlled method. A proposed method of "dumping properly" all 
particles emerging from the stripping foil( except the protons) will be discussed. 

4.1.2 Layout of the H Injection Region 

A schematic diagram of the Injection region is shown in Figure 14.11 In this diagram the 
following items are shown: 
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Figure 4.1: Schematic diagram of the H- Injection Region. 



a. Three of the main AGS combined-function main magnets (B19,B20, and C01). 

b. The stripping foil which is shown inside the magnet C01. 

c. Part of the closed orbit of the circulating proton-beam. The closed orbit is locally displaced 

"bumped" to the outside, and at the location of the stripping foil. The reason for "locally 
bumping" the closed orbit is partly described in Section 14. 1. li d) . e) above. Section 14.1.31 
below describes the method to generate the locally bumped orbit. At the end of the H~ 
injection time interval, the "locally-bumped" orbit of the circulating beam will return 
to a new orbit which is not bumped (labeled as "Circulating p-beam" in Figure 14.11 ) 

d. The trajectory of the injected H~ beam. The direction of the injected H~ is defined by the 

coordinate point (x,x') shown in the diagram. The coordinate x is the distance from the 
center point of the straight section SS_B19 and x' is the angle between the direction of 
the H~ and the SS_B19. At this distance x, the influence of the fringe field of the AGS 
main magnets on the H~is negligible. As the injected beam approaches the injection 
point (stripping foil shown in Figure 14.1(1 the fringe field and main field of the AGS main 
magnets B19 and C01 become significant and are taken into account in the calculations 
that determine: (1) the location of the stripping foil inside the magnet C01; (2) the 
beam parameters at the injection point (see Section 14. 1.41 below). 

e. The trajectories of the H~ beam that are not stripped by the stripping foil ( the partially 

stripped H~ beam and the electrons emanating from the stripping foil) must all be 
directed into a "dump" downstream of the stripping foil. 



4.1.3 The "Local Beam-Bump" at the Injection Region 

In order to make the circulating proton beam collinear with the H~ injected beam at the 

injection point (stripping foil), circulating proton beam is "locally bumped" by using two 
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horizontal "bump- magnets". The first "bump-magnet" is located at ~-90° phase advance, 
upstream of the stripping foil (location of the first "bump- magnet" is at the straight section 
SS_B12) , and the second "bump-magnet" is located at ~+90° phase advance downstream of 
the stripping foil (location of the second "bump- magnet" is at the straight section SS_C06). 
Figure 14.21 plots the displacement (y-axis) of the closed orbit at the center of the straight 
sections (x-axis) of the AGS for the two cases: 

a. When the "bumped-magnets" are turned off, and the orbit (red line) is not "locally 

bumped" ; 

b. when the "bumped-magnets" are turned on , and the orbit (blue line) is "locally bumped". 

The strength of the "bumped-magnets" whose location is systematically shown in Figure 
14.21 is ~3 mrad each. The horizontal and vertical tunes (Q x ,Qy) and the lengths of the 
closed orbits for the two cases is shown in the Figure 14.21 



at the midle of the AGS Straight Sections 
Backlog-windings; at B12 & C06 magnets 




R [nil 



Figure 4.2: The transverse displacement of the AGS closed orbit X co ^ at the middle of the 
straight sections of AGS. The red and the blue lines correspond to the X co d without bumps 
and with bumps respectively. 

Figure 14 .31 shows the displacement (y-axis) of a section of the "locally bumped" closed orbit 
(black line) at the location of the C01 main magnet of the AGS. The maximum displacement 
(~56 mm) of the closed orbit occurs at the center (~38 inches) of the 75 inches long C01 
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magnet. This displaced orbit has been generated using the two magnet mentioned earlier 
which are located at ±90° phase advance from the C01 magnet (this section). Each of the other 
lines are trajectories of the H~ injected beam corresponding to a different (x,x') direction of 
the H~ injected beam (see Figure 14.11 Figure 14.31 and Section 14.1.21 (d) for definition of x,x'). 

The values of (x x') that were used to generate the trajectories shown in Figure 14.31 are 
shown also on the same Figure 14.31 Each of the H~ trajectories is tangent to the "bumped- 
closed orbit" (black line) and this has been achieved by adjusting the initial conditions (x,x') 
of the H~ injected beam. 
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Figure 4.3: The black line corresponds to the "bumped closed orbit" X co d in the AGS magnet 
COl. This orbit has been generated using two magnet located at ±90° phase advance from the 
C01 magnet (see text). The COl magnet is 75 inches long. The rest of the lines are various 
trajectories of the H~ injected beam. Each of these trajectories corresponds to a different 
initial condition (x,x') of the injected H~ beam. 

From all the possible H~ trajectories shown in Figure I4.3[ the H~ trajectory that corre- 
sponds to the (x,x') = (250 mm, -56.5 mrad) has been selected as the "injection trajectory". 
Both, the closed orbit calculations for the derivation of the X co d as a function of distance, 
and the various H~ trajectories shown in Figure I4.3I have been calculated using the code 
AGS_BATE|TH|. This code is a modification of the code BEAM. The code uses the experi- 
mentally measured fields at the median plane of the AGS magnets to ray trace the charged 
particles. The code also provides transport matrices at any location along a trajectory and 
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beam parameters at any location along a closed orbit. Detailed description of the code and 
additional references are given in reference ^H]- in addition to the "local-bump" of the closed 
orbit discussed above, a set of two fast time varying bumps (horizontal and vertical) of the 
closed orbit my also be superimposed to help optimize the "beam painting" of the AGS phase 
space. 

4.1.4 Beam Parameters at the H Injection Point 

In order to minimize the emittance blow-up of the injected H~ beam, the beam parameters 
of the injected H~ beam must match those of the circulating beam at the "injection point", 
which is chosen to be the location of the stripping foil. 

The beam parameters of the circulating beam in AGS at the "injection point" have been 
calculated at the injection beam energy of 1.2 GeV, using the computer code AGS_BATE|15j 
and these beam parameters appear in Table 14.11 

The proximity of the H~ injection line to the AGS main magnet B20, shown in Figure 14 A\ 
suggests that the fringe field of the B20 magnet will affect the beam ellipse of the H~ injected 
beam. Therefore, the fringe field of the B20 magnet as well as that of the C01 magnet (see 
Figure I4.1|) are included in the calculations of the first order transfer matrices (R#,Ry) which 
define the beam transport between the "Matching point" (the point (x,x') shown in Figure 
14. II is defined in this note as the "matching point") and the "Injection point". 

Using the code AGSJ3ATE \Jj} we have calculated the first order horizontal and vertical 
transport matrices (R#,Ry) between the "Matching point" and "Injection point". These ma- 
trices appear below and correspond to the central trajectory with coordinate at the "matching 
point" (x,x') = (250 mm,-56.5 mrad). 

0.7015 6.058 
H ~ -0.0292 1.173 

0.7015 6.058 
v ~ 0.00497 0.796 

In order simplify the matching of the beam parameters of H~ injected beam to those of 
the circulating beam we have calculated the beam parameters at the "Matching point" (see 
Table 14. lj) by using both, the values of the beam at the "Injection Point", and the values of 
the transfer matrices R# and Ry above. 



Table 4.1: The beam parameters at the "Injection Point" and "Matching Point" (see text). 



Location 




a x 


rj x {m 




P v [m] 


a y 


Injection Point 


28.0 


-1.00 


-1.25 


-0.05 


8.00 


0.43 


Matching Point 


26.6 


-1.32 


-1.00 


-0.10 


16.10 


1.75 



4.1.5 H , H and Electron Beam Dump 

The H~ injected beam will not be stripped totally of its electrons by the stripping foil which 

is located at the injection point. Therefore the emerging beam from the foil will consists of: 
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Figure 4.4: Cross section of the B20 combined-function main magnet of the AGS, looking 
downstream. The oval line between the pole pieces is the vacuum chamber of the circulating 
proton beam. The circle shown between the top and bottom coils (red regions) of the magnet 
is an approximate location of the vacuum chamber of the H~ injected beam. 
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a. protons (H + ) which constitutes the circulating beam to be accelerated to 28 GeV. 

b. H~ ions that did not loose any electrons passing through the foil. 

c. H° neutral particles that lose only one electron. 

d. electrons (e~). 

Except for the fully stripped H~ injected beam, which is the circulating proton beam, the 
rest of the particles (H~, H° , and the electrons) must be "dumped" into a material which 
can absorb the deposited energy. The different rigidity of these three particles species, may 
require a specially designed "beam dump" . Similar problems of how to dump the partially 
stripped H~beam have been encountered in the 200 MeV H~injection in AGS[H], and the 
200 MeV H~injection in the Booster ^Hj- A study for the design of the beam dump of the 
partially stripped 1.2 GeV H~beam is required. 

4.2 Extraction from AGS at 28.0 GeV 

At present the AGS synchrotron is equipped with a Fast Beam Extraction (FEB) system 
which is used to extract the circulating bunched beam in AGS to the AtR beam transfer line. 
The FEB system will be used for the extraction of the 28 GeV high intensity proton beam 
which will irradiate the "neutrino production" target. A section (the U-line) of the AtR beam 
transfer line which is used to transfer the extracted beam from AGS into the Relativistic 
Heavy Ion Collider (RHIC) will also be used to transport the 28 GeV proton beam to the 
"neutrino production" target. Detailed description of the FEB system has been presented in 
previous technical notes therefore in this report we will only outline the main features 

of the (FEB) system. 

4.2.1 Layout of the Extraction Region 

The main steps of the procedure for fast beam extraction are outlined below: 

a. Beam acceleration from injection energy to 28 GeV. 

b. Generate two "local beam bumps", one (BLWG09 see Figure 14. 5J) which displaces the 

circulating beam into the C-type Fast kicker magnet (FK) which will "kick" the beam 
out of the AGS, and the other (BLWH11 see Figure 14.51 ) which displaces the circulating 
beam close to the H10 septum magnet (SM) which is located downstream from the G10 
kicker at ~270° phase advance. It is advantageous to use a C-type Fast kicker magnet 
which is displaced from the ideal orbit, instead of a full aperture kicker 

c. The circulating beam bunches pass through the aperture of the G10 Fast kicker, and are 

kicked by the G10 Fast kicker which displaces the beam bunches on the other side of 
the H10 septum, into the U-line (see Figure 14.51 ) which will be used as the first section 
of the beam transport line to the "Neutrino production" target. A schematic diagram 
of the FEB system in AGS is shown in Figure 14.51 In this figure shown are: 
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(i) the BLW_G09 and BLW_H11 "local-bumps"; 

(ii) the Fast Extraction kicker (FK); 

(iii) the Beam Extraction Septum Magnet (SM); 

(iv) and the beginning of the U-line. 

Detailed description of the FEB system is given in references |H)],and |17j . 

4.2.2 Extraction Devices 

In this section we describe briefly the function of the devices which are used for the FEB in 
the time order the devices are involved in the fast beam extraction. 

1. The two 3A/2 "local beam bumps" BLWG09 and BLWH11 (see Figure Each of 
the bumps is generated by back-leg (BLW) windings on the return iron of specified main 
magnets of the AGS. 

a. The BLWG09 "beam-bump" which displaces the closed orbit inside the aperture of 

the G10 fast kicker that is displaced radially from the ideal orbit on the outside 
of the AGS ring, because it has acceptance large enough to accept the size of the 
28 GeV beam but not large enough to accept the size of the injected beam. A 
proposed "full-aperture" kicker which can accept the size of the injected beam, 
does not need to be displaced radially nor does the beam need to be bumped by 
the BLWG09 into the kicker. The large gap of a "full aperture" kicker, however, 
will increase the power supply requirements of the fast kicker. 

b. The BLWH11 "beam-bump" which displaces the closed orbit of the circulating beam 

near the septum of the H10 extraction magnet (SM). This displacement of the 
beam closed to the septum magnet reduces the strength requirements of the kicker 
to displace the beam into the other side of the septum and into the aperture of the 
magnet. 

2. The G10 fast extraction kicker (FK Figure 14.51 ) which kicks each circulating beam- 
bunch into the aperture of the septum magnet. 

3. The H10 extraction septum magnet (SM Figure 14.51 ) which extracts the beam bunch 
from the AGS into the U-line. 

4.2.3 Beam Parameters at the Extraction Point 

The extracted beam to be transported to the "neutrino production" target will use part of 
the AtR transport line, the U-line which will be extended to the production target. The beam 
transport line is discussed in more detail in one of the chapters of this report. In this section 
we only report the beam parameters of the extracted beam (Table 14.21 at the extraction point 
of the transport line. The knowledge of the beam parameters of the extracted beam at the 
beginning of the transport line are required for the design of the transport line. Details of 
the method which is used to calculate the beam parameters and the calculation of the beam 
parameters are given in references ^H] and |19j . 
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Figure 4.5: Schematic diagram of the AGS RING, Fast-Beam-Extraction (FEB) system. 
Shown are: a) the BLW_G09 and BLW_H11 "local-bumps" b) the Fast Extraction kicker 
(FK) c) the Beam Extraction Septum Magnet (SM) d) and the beginning of the U-line. 



Table 4.2: The beam parameters at the "Extraction Point" of the beam transport line to the 
"neutrino production" target (see text). The dispersion and angular dispersion (rjxrfx) were 
not measured. 



Location 


&M 


a x 


%H 


r}\ 


P y [m] 


Oiy 


Extraction Point (Theory) 


37.5 


-4.1 


-1.5 


-0.13 


6.5 


0.85 


Extraction Point (Measr.) 


36.7 


-4.5 






7.7 


1.2 
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5 Beam Transport to Target 

The extracted beam will enter into the existing U-line of the AGS. For aiming the neutrino 
beam to the Homestake site in South Dakota, the proton beam has to bend 68 degrees and 4 
seconds horizontally and 11.26 degrees downwards. Figure 15.21 shows the proposed layout of 
the beam transport to the neutrino target at the Homestake site. 
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Figure 5.1: Elevation view of the neutrino beam line to Homestake, South Dakota. 



To aim 11.26 degrees downwards, the proton beam has to climb up a high hill to the target 
and decay channel. The hill arrangement keeps the target and hadronic decay channel well 
above the water table of Long Island. Figure 15.11 shows the sketch of the hill. Figure 15.31 
shows the cross section of the hill. 

The beam is bent vertically upward by 11 degrees with help of one cell with phase advance 
of 90 degrees. In the vertical incline plane, beam is then bent 68 degrees and 4 minutes hori- 
zontally in four cells with a phase advance of 90 degrees. Finally, the beam is bent downwards 
by 11.26 degrees; aiming it towards Homestake, South Dakota. Before the final bend, the line 
uses one doublet to focus the beam with 2 mm radius (rms) on target. This arrangement 
results in an achromatic system in both horizontal and vertical planes. The transport system 
uses a 4Q16 existing quadrupole and dipole for the vertical bend at 1 degree/meter and for 
the horizontal bend at 0.5 degree/meter. Figure 15.41 shows the (3 and rj functions for the beam 
transport to the target from the U-line (ULTBT). 
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Near Detector 

<]=] NORTH 

Figure 5.2: A plan view of the neutrino beam facility at the AGS. The 1 MW proton beam 
will be taken from the existing fast extraction line (the U-line) and continued on the up-slope 
of the hill. The vertical and horizontal bending magnets are separate in this plan. The beam 
will be bent downwards at the top of the hill to aim it towards Homestake at 11 degrees. 
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Figure 5.3: Hill cross section at target area. 



U-LINE TO TARGET TRANSFER LINE 




S[m] 

Time: Thu Apr 1 1 4:1 9:47 2003 Last file modify time: Thu Apr 1 1 4:07:21 2003 

Figure 5.4: (3 and 77 functions along the ULTBT. 
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6 Target/Horn Conceptual Design 
6.1 Background and Overall Issues 

To achieve the 1 MW upgrade option of the proton driver at BNL serious consideration must 
be given to both the target material and selection and horn configuration. A solid target is a 
viable choice for a 1 MW beam. Low and high Z materials have been investigated both in 
terms of the material endurance as well as the feasibility of target /horn configuration. Results 
of the parametric studies on material choices regarding pion production are shown in 16.11 
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Figure 6.1: Pion production as function of proton momentum for different materials. 

A graphite-based carbon-carbon composite is selected as a target material both for its low-Z 
and thermo-mechanical properties. Various aluminum-based materials are being considered for 
horn conductors, such as 3000-series and 6061-T6 aluminum, that will combine low resistivity, 
high strength and resistance to corrosion and micro-cracking. According to the baseline design 
the following beam parameters are guiding the material choice, horn and target design, cooling 
schemes, etc. 

In the process of designing the horn/target system the following key elements were closely 
addressed: 

• Heat generation and removal from the target/horn system 

• Target thermo-mechanical response from intercepting energetic, high intensity protons 

• Irradiation and corrosion effects on materials 

• Horn/target integration issues 

• Horn mechanical response and long term integrity (irradiation, corrosion and thermal 
fatigue) 
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Table 6.1: Horn/Target design parameters. 



Proton Beam Energy 


28 GeV 


Protons per Pulse 


8.9xl0 13 


Average Beam Current 


35.7 /iA 


Repetition Rate 


2.5 Hz 


Pulse Length 


2.58 /isec 


Number of Bunches 


23 


Number Protons per Bunch 


3.87xl0 12 


AGS Circumference 


807.1 m 


Bunch Length 


40 ns 


Bunch Spacing 


60 ns 


Normalized Emittance-X 


100 7r mm-mrad 


Normalized Emittance-Y 


100 7r mm-mrad 


Longitudinal Emittance 


5.0 eV-sec 


Target Material 


carbon-carbon composite 


Target Diameter 


1.2 cm 


Target Length 


80 cm 


Horn Small Radius 


7 mm 


Beam Size (Radius) on Target 


2 mm (rms) 


Horn Smallest Radius 


6 mm 


Horn Large Radius 


61 mm 


Horn Inner Conductor Thickness 


2.5 mm 


Horn Minimum Thickness 


1 mm 


Horn Length 


217 cm 


Horn Peak Current 


250 kA 


Current Repetition Rate 


2.5 Hz 


Power Supply Wave Form 


Sinusoidal, Base Width 1.20 ms 
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• Beam windows integrated in the system to (a) separate the vacuum space in the transfer 
line from the final beam run to the target and (b) to maintain the coolant around the 
target in a close-system loop. 

6.2 Description of the Integrated System 

Figure 16.21 is a conceptual description of the target and horn integrated system being con- 
sidered in this study. The 1.2 cm diameter carbon-carbon composite target is fully inserted 
into the inner horn conductor allowing a 1 mm annular gap between the target and horn 
surfaces for the helium coolant to flow and remove heat. Shown in the front of the target is 
a carbon-carbon beam "collimator" or baffle that has dual role. Specifically, it provides the 
upstream target support and includes the channels that will force the coolant flow into the 
annular space. Further, it plays the role of the beam diffuser in the event the proton beam 
strays off beam axis and thus protects the horn. The size of this baffle needs to be optimized 
to play that role effectively. In addition to the above two functions, the front end of the target 
will be maintained at a low temperature which will help in removing heat deposited on the 
target by conducting into the mass of the baffle. 

It is envisioned that the baffle/target arrangement will be a monolithic structure for system 
optimization. At the downstream end of the target, a special design allows for the forced 
coolant to leave the horn's neck-down section and also provides a simple support for the 
target in the current design. The option of maintaining the target as a cantilever beam (only 
supported from the upstream end is being considered. 

The horn, made out of an aluminum alloy, has a diameter in its narrowest section of 1.4 
cm while the thickness over that section is 2.5 mm. The thickness of the inner conductor 
reduces to ~ 1 mm downstream of the neck-down section. As shown, the current supply 
and return take place at the downstream end of the horn. The baseline design calls for a 
250 kA peak current to be achieved with a half-sine wave shape that has a base of 20 /xs 
with a repetition rate of 2.5 Hz. The magnetic pressures and joule heating generated in 
the conductor are the driving forces in the optimization of the horn conductor. While heat 
generated in the narrowest section of the horn will be partly removed by the fluid flowing in 
the annular space, heat generated in the skin of the current surface will be being removed 
by spraying a coolant through a set of optimally positioned jets. Two options are being 
considered, namely, the spraying of water and cold helium gas. The schematic of Figure 16.21 
depicts the arrangement with water-cooling. The choice of cooling, as discussed later in detail, 
is significant in maximizing the useful life of the aluminum material in that environment. 

Under consideration, and shown in Figure 16.21 is a downstream thin window whose role 
is to hold the target coolant in a closed system. This way the coolant is collected, cooled 
and returned to the target upstream to be re-ejected into the annular space. The key issue 
with such window is the fact that it will see a significant portion of the incoming beam power 
(~ 80%) and will be subjected to high thermo-mechanical stress conditions. Further, the 
presence of additional material in the flight path of pions generated and focused by the horn 
is a further impediment. However, since the only role of such window is the prevention of 
coolant to escape from the closed envelope, the materials considered are low-Z (for minimal 
interaction with secondary particles or generation from intercepting the beam protons) and 
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minimal heating, which combined with the small thickness will minimize the problem further. 



l_out 




Carbon-Carbon 
Target Rod 



Figure 6.2: Conceptual arrangement of horn and target. 



6.3 Evaluation of Heat Removal Principles 

As part of the design process, a number of cooling schemes have been considered for the 
integrated system of horn and target. Utilizing past experience on horn/target system per- 
formance we arrived at the current baseline scheme which calls for a de-coupling of the target 
from the horn with regards to their cooling and the use of helium gas as the fluid to remove 
the heat from the target. 

Specifically, as shown in Figure 16.21 cold helium is forced past the surface of the CC 
target at sufficient velocity to ensure that the heat generated per beam pulse is removed 
before the next pulse arrives. The helium flowing in the annular space is also responsible 
for removing some of the heat generated in the horn inner conductor that surrounds the 
target. This heat, as will be discussed in detail, originates in three separate mechanisms, 
namely, joule heating from the flowing current, energy deposition due to the proton beam 
interaction with the target, and heat by radiation from the surface of the CC target. Further, 
heat generated in horn conductors (inner & outer) due to current and beam is removed with 
spraying of cooling water over the inner horn surface. Since different mechanisms, namely 
heat convection, conduction and radiation between target and horn, will be at work, the heat 
balance of the overall system as it reaches an operating temperature will be addressed by the 
use of a sophisticated finite element model which is described in a later section. As a first step, 
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however, and in order to assess the heat removal capacity of surfaces in the system, a set of 
basic calculations are performed that will indicate whether the adopted schemes are feasible. 
Furthermore, the results of this preliminary analysis will serve as a benchmark and check for 
the finite element simulations. The following analysis addresses the joule heating in the horn, 
the radiation heat exchange between the target and the horn conductor, the heat removal 
capacity of the flowing He in the annular space and, the heat removal capacity/requirement 
of the water spray inside the horn. 

6.3.1 Joule Heating Estimate in the Horn 

Given the current pulse structure shown in Figure 16.31 with an effective frequency of 0.025 
MHz, the current is expected to flow over a skin depth of the inner surface of the inner 
conductor. Key to the evaluation of joule heating is the estimation of the skin depth for the 
parameters of this study. 

The skin depth 5 for a horn conductor made out of 3000-series aluminum, for example, 
with resistivity p = 4.2 mohm-cm, is calculated based on the following relations: 

6 = (em/f^k, (6.i) 

where, f = 0.025 MHz, k 1= [p/p c ] 1/2 and p c= 1.724 mhoms -cm is the resistivity of copper. 
From these relations one obtains, 

S Al = 0.06525cm . (6.2) 




time 



Figure 6.3: Current pulse structure. 

By focusing attention on the neck-down section of the inner conductor where conditions are 
most severe, the joule heating generated by the current pulse and the anticipated temperature 
rise is estimated. The current density is a function of the penetration into the conductor 
according to the relation: 

J{z) = J e~ z/5 (6.3) 
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such that, 

I{z) = I J(z)dA, (6.4) 



A 



where, J D = 689 kA/cm 2 is the current density at the conductor surface and z is the 
penetration and dA = 2tt (R, out — z) dz, where H out = 0.95 cm. The joule heating generated 
by the pulse per cm-length of the inner conductor is 



20fis 

J H n _ 

A 

where, J(z,t) = J(z) sm(irt/20fis). By carrying out the integration above one obtains, 



fp,„«-™= J jAz,t)pdAdt. (6.5) 



JH/ (pulse-cm) = 3.88 Joules. (6.6) 
The total power generated in the neck-down section of the horn alone is estimated to be: 

P neck = 2.5(Hz) ■ 3.88 J/(cm-pulse) • 86 cm = 835 Watts. (6.7) 

The peak temperature rise in the conductor, experienced at the inner surface where the 
current density is at maximum during each current pulse, is estimated to be 



AT = 7.4°C. 
6.3.2 Estimate of Radiation Heat Transfer 

Given that the CC target will be operating at much higher temperatures than those of the 
horn (~ 840 °C) some heat will radiate from the target surface to the surface of the horn over 
the neck-down section. Assuming that the surface temperature of the aluminum is maintained 
at ~ 90 °C, then the radiating heat flux can be estimated from the relation: 

qcc->Ai = cr [Tec ~ T Al] /g 
A C c _L_ _i_ Aon ( J_ _ A' 

£cc Aai \£ai J 

where, a = 5.669 x 10~ 8 W/m 2 - K 4 , e cc = 0.98, e A i= 0.19 

Based on the above parameters, the heat flux from the target surface is approximately q/A 
= 11916 W/cm 2 and the total heat transfer from the target to the horn inner conductor q— 
11916. A target = 360 W. As seen the amount of heat exchanged by heat transfer is quite 
low. However, if the surface of the horn seeing the target is treated in order to increase its 
emissivity to Eai= 0.19, then the heat transferred to the horn amounts to ~ 1.363 kW. 



59 



April 15, 2003 



AGS Super Neutrino Beam Facility 



Heat Removal from Horn Neck-down 



6.3.3 Target Cooling Calculations 

Given that the amount of heat transferred by radiation alone from the target to the horn is 
quite small, the primary mechanism is heat removal by forced convection of cold He in the 
annular space between the target and the horn. The key elements in assessing how feasible such 
a mechanism is are: the operating temperature of the target; the temperature of the cooling 
fluid; and the heat generation rate in the volume of the target that needs to be removed. 
These working parameters will establish the velocity of the coolant necessary to remove the 
generated heat. Using as a baseline, the energy deposition in the volume of the target for a 2 
mm RMS beam spot, (which amounts to 7.1 kJ/pulse or 17.75 kW) as an acceptable average 
operating temperature of surface of the target T s = 800°C and a temperature of the helium 
of ~ 144 K, the flow velocity is estimated as follows, 

Heat flux from target surface = q/A tar get = hf(T a — X# e ). (6.9) 

One can estimate the necessary film coefficient hf from the above relation to be 618 W/m 2 -C. 
By relating the film coefficient to the Nusselt number Nuo e 

Nu De = h f D e /k = 13.326, (6.10) 

where, k = 0.14 W/m-K and D e = 0.014-0.012 = 0.002 m (equivalent annular diameter). By 
relating the Nusselt number to the Reynold's number N# e through the relation 

Nu De = 0.023(iV fe ) a8 (P r ) a3 , (6.11) 

where P r is the Prandtl number, the Reynolds number is estimated to be ~3400. Reynolds 
number on the other hand is related to the velocity of flow Uy through, 

N Re = U s D e /u. (6.12) 

where v is the dynamic viscosity of He iy = 37.1 x 1 0~ 6 m 2 /s). From the above relations it 
is estimated that a flow of He with velocity U/ = 65 m/s is required. 

These calculations, however, do not account for the heat moved away from the target by 
radiation and conduction. The later takes place near the front of the CC target by maintaining 
the block of the baffle (or collimator) that holds the target in place at a low temperature. 

6.4 Heat Removal from Horn Neck-Down Through Forced Convec- 
tion 

In this section estimates are made of the required convective heat transfer capacity of He 
in the annular space to enable the removal of all the heat generated in the volume of the 
86cm-long neck-down section. The best estimates, thus far, of the heat deposited over that 
volume follows: 

Joule Heat = 1.335 kW 

Heat from beam interaction with target = 36% of 17.75 kW = 6.39 kW 
Radiation from target = 1.363 kW 
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TOTAL = 9.088 kW. 

The high value of energy deposited in that section alone from beam-target interaction needs 
further verification. Independent analysis is currently being performed to verify this energy 
deposition. The available area of the horn surface with radius Rj = 0.7 cm is Aj = 0.0352 m 2 . 
This leads to a surface flux of 

q/Ai = 258, 181W/m 2 . (6.13) 

With forced helium temperature of 144 K and surface temperature of the horn ~ 90° C, 
the required film coefficient is hf = 1179 W/m 2 -C. Following the procedures of the previous 
section that relate the convective film coefficient with the Nusselt and Reynolds numbers, one 
arrives at a flow velocity required to remove the heat in the neck-down section equal to U# e 
= 135 m/s. 

6.4.1 Heat Removal from Horn Inner Surface by Coolant Spray 



Option I - RAW Spray 

The baseline method utilized to remove the generated heat over the entire inner surface (cur- 
rent side) as well as the energy deposited throughout the horn from beam-target interaction 
is by spraying water through special nozzles that penetrate the outer conductor. While this 
approach has been used extensively to cool horns in the past (and in many instances to cool 
the target as well through conduction across the horn inner conductor), there are some issues 
associated with such approach, namely, the potential acceleration of corrosion of the aluminum 
surface coming in contact with water and the fatigue failure through propagation of surface 
cracks plus the fact that water within the horn will be in the path of pions. The attempt 
to maintain the aluminum surface temperature, when in contact with water, below 100° C 
(based on experience from the reactor industry) in an attempt to extend the life of the horn, 
results in larger quantities of water being sprayed against the conductor inner surface. 

To estimate the required spraying capacity, attention is again focused on the inner con- 
ductor which will experience higher joule heating and gamma-ray heating. Assuming that the 
spraying jets are positioned in such a way that the entire inner surface experiences forced flow 
sweeping by (resembling a cylinder in cross- flow), the following forced convection relations 
apply, 



Nu D = 0.3 



0.62Re% 2 Pr 1/3 



1 + 



282000' 



4/5 



(6.14) 



[1 + (0.4/ Pr) 2 / 3 ] 1/4 

where, Nup = hD/k is the Nusselt number and h is the effective film coefficient. 

Assuming that the water jets force a flow with a free velocity of U = 2.5 m/s at a tem- 
perature of T wa ter = 20° C and that the conductor surface is maintained at a temperature of 
Twaii = 80° C, the following fluid properties apply: Reynolds number Re^ = U D/v = 87963, 
Pr = 3.6, v = 0.0054 cm 2 /sec and k = 0.645 W/m-K. 
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Substitution into the above relation leads to, 

Nud = hD/k = 367. (6.15) 

And to an effective film transfer coefficient of: h = 12485 W/m 2 -K . 

By relating the heat flux from the conductor surface to the convective heat transfer, 

q/A = h{T wall - T water )W/m 2 , (6.16) 

the heat transferred through the surface area A is q = 38.5 kW. Clearly this cooling scenario 
can satisfy the needs for removing the heat from the inner surface of the conductor even for 
operations that will increase the heat deposition in the conductor. 

OPTION II - Helium Spray 

In an effort to eliminate the RAW altogether, an attempt is made to remove excess heat by 
spraying He onto the surfaces instead of water. To estimate the capacity of such system, we 
again focus on the neck-down section that represents a horizontal cylinder of radius R G = 9.5 
mm and length L = 86 cm. The relations that apply are those of gases flowing normal to 
a single cylinder. The Reynolds number for this configuration is N^ e = U m D D /V while the 
Nusselt number is related to it by the empirical relation 

Nu Do = C(N Re ) m . (6.17) 

Where C and m are parameters that are based on the range of the Reynolds number. For 
a free velocity of He (after leaving the nozzle and before impinging on the surface) of 30 m/s 
and temperature T# e = 144 K, while maintaining the aluminum surface temperature to 90° 
C, one can estimate N^ e = 1536 and, 

Nu Do = 54.9 = h f D /k. (6.18) 

Which leads to h f = 270 W/m 2 -C. 

Using the relation: q/A G = hj (T D - T# e ) where A D = 0.0513 m 2 , the heat transfer capacity 
of the inner surface is approximately 3000 W. Any increase of the velocity of He coming out of 
the nozzle or a potential reduction in the heating of the horn from beam-target interaction, will 
make the cooling scheme that depends entirely on He sufficient to remove the heat deposited 
on the horn. Further experimental studies are needed to assess the potential of this option. 

6.5 Target Consideration 

Material Selection — Issues 

The baseline material for the superbeam experiment is a special Carbon-Carbon composite. 

It is a 3-D woven material that exhibits extremely low thermal expansion for temperatures up 
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to 1000 °C. For higher temperatures it behaves like typical graphite. Its thermal expansion 
behavior is significant in the sense that the thermo-elastic stresses induced by intercepting the 
beam will be quite small, thus extending the life of the target. The Table 16.21 below lists the 
variation of expansion as a function of temperature. 

Table 6.2: Thermal expansion data for a 3-D fine weave Carbon-Carbon composite. 



Temp. 


% Elongation 


23 °C 


0% 


200 °C 


-0.023% 


400° C 


-0.028% 


600° C 


-0.020% 


800° C 


0% 


1000° C 


0.040% 


1200° C 


0.084% 


1600° C 


0.190% 


2000° C 


0.310% 


2300° C 


0.405% 



In addition to the low thermal expansion, the material is stronger than typical graphite, 
reaching compressive strengths > 120 MPa. 

Experience on this material has been acquired as a result of the BNL E951 experiment 
and its use in the SNS foil scrapers. Some of the findings, relevant to this study, are discussed 
in a later section. 

The long term behavior of this special material in an irradiation environment is not known. 
There was an initiative, as part of the E951 experiment, to irradiate this material to dpa levels 
that could be comparable to running times that are approaching anticipated lifetimes for this 
experiment. While this particular material was not irradiated, it is anticipated, as part of 
the R&D effort, to be tested in the future to enable the assessment of its integrity after the 
accumulation of several dpa. 

Beam Energy Deposition and Heat Removal 



In the current option, the target is an 80 cm long cylindrical rod with two competing 
diameter sizes (6 mm and 12 mm). The 6 mm diameter target is chosen to intercept the 100- 
TP 1 mm rms proton beam and the 12 mm for the 2 mm rms beam. There are advantages and 
dis-advantages in both choices. When intercepting the 1 mm beam, the total energy deposition 
in the target per pulse amounts to about 5.1 kJ. This corresponds to a peak temperature rise 
of about 1000 °C in the center of the target. With the larger beam spot (2 mm rms) the total 
energy deposited and converted to heat in the target is 7.3 kJ per pulse with peak temperature 
rise of about 270 °C. Heat will be removed from the target through forced convection of helium 
past the outside surface. While the 2 mm beam/12 mm target diameter option appears to be 
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the optimum, the fact that the heat within the target will have to travel a longer distance to 
reach the surface and be removed, may make the 1 mm beam option attractive. 

Engineering calculations and detailed finite element analysis were employed to qualify the 
two options while taking into consideration the special target/horn configuration. Figure 16.21 
depicts the current thinking of system layout. Specifically, the target rod is placed within 
the neck-down section of the inner horn conductor forming an annular space with the inner 
surface of the conductor. The long target is supported at both ends. In the front, a specially 
machined carbon-carbon block holds the target while at the same time it allows for the beam 
to go through without intercepting it. In addition the this special block allows for special 
nozzles that will force helium into the annulus for target cooling and it provides a conduction 
path in the vicinity of the target section that is heated the most. Further, the block will play 
a diffusing role in the event of a beam straying from center. At the far end, the target includes 
four special fins that allow it to rest on the inner surface while allowing the helium to pass 
through. A special notch on the inner conductor will prevent the target from moving further 
into the horn. 

Calculations of heat removal of independent heat transfer mechanisms while using the basic 
principles were performed and are presented in an earlier section. The engineering calculations 
for removing the deposited energy revealed the following: Since the target is to be decoupled 
from the horn, the forced helium (and the conduction near the front of the target) must remove 
all the heat that amounts to 18.25 kJ for the 2 mm rms beam/12 mm diameter target. The 
energy deposition varies along the beam and so the engineering calculations based on the total 
energy in the volume can only provide a rough estimate. Analytic calculations are expected 
to refine the engineering numbers. Specifically, to remove the 18.25 kJ from the surface of 
the target through forced convection, a helium temperature of ~5 °C and a target operating 
temperature of 840 °C were assumed. The heat transfer calculations (using detailed finite 
element analysis) reveal that helium flowing continuously in the annulus with velocity near 
40 m/s will, in conjunction with radiation to the horn surface and conduction near the front 
of the target, enable the removal of all the heat deposited in the target. 

The detailed finite element model incorporated both the target and the horn in one single 
system. The only connection between the horn and the target was provided by the radiation 
exchange of the target surface and the surface of the horn inner conductor. Further, the 
two surfaces shared similar values of convective heat transfer due to the flowing helium in 
the annulus. In order for the overall system to achieve a "steady-state" (or base operating 
temperatures) throughout a transient analysis with all the heat transfer mechanisms in place 
was performed. 

Shown in Figure 16.41 are the transient temperatures at the hottest point in the center of 
the CC target and at the surface that result from a 100 TP/28 GeV/ 2 mm RMS proton beam 
on target. As mentioned previously, the temperature rise per pulse in the target is of the order 
of 270° C. For a 1 mm RMS beam spot on a 6 mm diameter target the temperature rise per 
pulse is ~ 1000° C, as shown in Figure 16.51 The total energy deposited on the target, however, 
is only 5.1 kJ per pulse as compared to the 7.3 kJ for the 2 mm RMS beam. While the benefit 
in smaller energy deposition would be welcomed, the temperatures in the CC target will be 
in the regime of higher strains and the resulting stresses will exceed the material strength. 
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Transient temperatures in the CC target (T_cl — target center ; T_s — target surface) 
intercepting a 100 TP/28 GeV/2nun lUVIS proton beam. Target diameter = 1.2 cm 

150 Ot 1 1 1 1 1 1 1 1 1 1 




250 I 1 1 1 1 1 1 1 1 1 

□ 4 8 12 16 20 

2 6 10 14 10 
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Figure 6.4: Transient temperatures in the CC target subjected to a 100 TP/28 GeV/2 mm 
RMS beam. Coolant temperature in the annular space T# e = 5 °C. 
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Figure 6.5: Temperature rise per pulse in the 6 mm-diameter CC target from a 100 TP/28 
GeV/1 mm RMS beam. 
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Evaluation of the Option with the Target & Horn in Contact 



A scheme used in past experiments involves solid targets and horns and allows for the 
removal of heat from both the target and the horn by cooling the outer face of the horn 
conductor (current side) through massive water spray. In this arrangement, the solid target is 
housed within the narrowest inner conductor and through surface contact it passes the heat 
generated by the beam/target interaction onto the conductor. The heat travels by conduction 
to the opposite face (current side) where it is removed by the forced cooling. To assess its 
potential application in this experiment, the system was analyzed by considering a perfect 
contact between the two systems (CC target and aluminum horn) and the required heat 
removal capacity was estimated. 

In performing this analysis, the following key points are kept in mind, namely: 

• Energy deposition and thus heat generated in the target in this study are almost an 
order of magnitude higher than the previous exercises; 

• The repetition rate in this study is also higher; 

• The goal in this study is to hold the aluminum temperatures as low as possible for 
surfaces that are in contact with water. 

An additional concern is the fact that, at the operating temperatures of the target, it will 
be impossible to maintain perfect contact between the surfaces. This will inevitably lead to 
the formation of hot spots and the subsequent melting of the aluminum. 

Figure 16.61 depicts the two materials in contact as they were modeled in the finite element 
analysis) and the temperature distribution when the aluminum surface is being cooled with 
an effective film coefficient of 1000 W/m 2 -C. As seen in Figure 16.61 a large section of the 
aluminum conductor (only a segment is shown in the figure) operates at temperatures near the 
melting point of aluminum (660° C). A parametric study revealed that in order to maintain 
temperatures in the aluminum that are below 100 C would require a film coefficient of 3000 
W/m 2 -C. Such a heat transfer coefficient will be almost impossible to achieve in such setting. 
However, even if that was achievable, the issue of hot spot formation cannot be eliminated 
since the two materials expand differently. Thus this cooling method is not suitable here. 

Thermo-elastic target Stresses 



Figures 16.71 and 16.81 depict the dynamic stresses in the target intercepting the train of six 
micro-bunches that comprise each 100 TP pulse. The resulting equivalent stresses, estimated 
for a pulse which arrives when the target has achieved a "steady-state" operating temperature 
(~840° C) peak at ~ 38 MPa. Such stresses are well below the strength limits of ~ 120 MPa. 
However, as it will be discussed in the following section, experimental measurements have 
shown that the material exhibits some dynamic inertia during fast beam energy deposition 
that is expected to slightly increase the anticipated stress levels that are based on the thermal 
expansion properties which, in turn, are estimated with slow heating of the material. 
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Figure 6.6: Evaluation of cooling scenario with elimination of cooling in the annulus and heat 
removal by conduction through target/horn contact. 
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Figure 6.7: vonMises dynamic stresses generated in the 2.2 cm diameter carbon-carbon 
composite and their attenuation before the arrival of the subsequent beam micro-bunch train. 
Results depicted for a 6-microbunch pulse train totaling 100 TP. 
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As a result of beam optimization, the new parameters call for a beam spill of 23 40 ns-long 
micro-bunches spaced by 60 ns. While calculations for such case have not been performed, 
it is expected that the thermal shock stresses will experience some relaxation between micro- 
bunches thus driving the stresses even lower. 
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Figure 6.8: Radial and axial stresses in the 1.2 cm diameter carbon-carbon target intercepting 
a 6-microbunch train 28 GeV proton beam Totaling 100 TP. 



Assessment of Long-term Target Survival 



Given that little is known about the behavior of this target material when it is highly 
irradiated, it is envisioned that, as part of the R&D, representative target material will be 
exposed to levels of dpa anticipated to be seen during the lifetime of the target and the its 
long-term integrity be assessed. 

Verification of CC Mechanical Properties 



While the chosen operating temperatures (including the temperature rise) are still within 
the very low thermal expansion regime and the expected thermal stresses and strains are within 
the safety envelope, some questions remain regarding the response of the target material 
during fast deposition of heat. In the course of experiment E951, 1 cm diameter carbon- 
carbon targets (along with ATJ graphite) were exposed to the 24 GeV AGS beam and their 
response was recorded on the outer surface in the form of axial strains. Figure 16.91 below 
shows the target arrangement. Figure 16. 101 depicts strains recorded on the surface of the CC 
target along the axial dimension. While these strains are, by comparison, much smaller than 
those recorded on the ATJ graphite target with same beam parameters (Figure I6.11j) . the 
generated temperatures from the beam energy deposition were very small to cause these strain 
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amplitudes (based on the material property data sheet). It appears that the material responds 
dynamically during very rapid energy deposition. For the 2 mm rms beam the anticipated 
thermo-mechanical stresses are well below the material strength limits even with a potential 
increase in stresses due to inertia in the material (as shown in Figure I6.10j) . It is anticipated 
that further tests of intense proton beams on such targets will shed light in verifying the true 
thermo-mechanical response of this special material. 

Figure 16. 121 depicts the predictions in axial strain on the AT J graphite target at a location 
similar to the location of the two fiber-optic strain gauges of Figure 16.111 The predictions were 
made using similar transient finite element methodologies as the one utilized in the analysis 
of the CC target and horn. It is clear that these computational techniques can provide an 
excellent estimate of the actual system response. 




Figure 6.9: Carbon-Carbon composite and ATJ graphite target arrangement for the E951 
experiment. 



6.6 Horn Design Issues 

Figure 16.21 depicts the first horn made of aluminum and with a thickness in the neck-down 
section of the conductor is 2.5 mm. The radius of the inner conductor for the following 
calculations was assumed to be 7 mm. The issues associated with the horn are the following: 

• Joule heating in the inner conductor 

• Energy deposited on the horn due to beam-target interaction 

• Heat received from target by radiation exchange 



Heat removal for the 2.5 Hz operation 
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Figure 6.10: Axial surface strain measurements in the CC target. 




Figure 6.11: Axial strains measured in the ATJ graphite target during E951 experiment. 
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□ 80 160 240 320 400 



Figure 6.12: Predicted dynamic strains at approximately same locations as strain gauges. 

• Thermal stresses (including highly transient loads) 

• Magnetic field pressures 

• Vibration of thin horn conductors (from pulsed load and possibly from He coolant) 

• Differential expansion between inner and outer conductor 

The most serious issue is the amount of heat generated and deposited in a small volume 
of the inner conductor where current is flowing mostly in the skin depth. For the engineering 
calculations 240 kA will flow into the horn during each pulse in a half-sine form with duration 
of 20 /is. 

Skin Depth — Effective Resistance — Joule Heating Calculations 

In a previous section the parameters governing the current flow in the horn and the re- 
sulting joule heat were estimated. The primary goal, given the structure of the current pulse 
is to estimate the skin depth and subsequently the joule heating throughout the horn. For 
a horn made of aluminum 3000-series and the 20 ps half-sine 240 kA peak current pulse the 
following relation lead to the estimate of skin depth and joule heating: 

• Resistivity p of Al-3000 = 4.2 mohm - cm 

• Skin depth 5 = [6.61/f 1 / 2 ] k x 

. k! = [p/p^ 2 



p c = 1.724 mohms-cm 
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Utilizing only the dominant frequency (f = 0.025 MHz) the skin depth has been estimated 
to be: 

S al =0.06525 cm. 

Using the parameters above (and a 240 kA current in the horn), calculations were made 
with the simplified assumption that all the current flows uniformly over one skin depth. The 
purpose of this simplification was to generate joule heat and temperature rise estimates that 
can be used to validate and benchmark the finite element simulation model being employed 
to estimate the heat generation and transfer in the integrated system. By assuming that the 
current only flows over one skin depth, an effective resistance can be estimated, 

R eff = p/A s = 4.2e - 6 (Ohm-cm) /(2tt • 0.95 • 0.06525) = 0.107836 ■ 10" 4 Ohm/cm. (6.19) 
For this simplified case the joule heating per current pulse can be estimated from 

JH pu i se - cm = J I 2 R e ffdt, (6.20) 
o 

where, I— 250 kA sm(7rt/20ps). From the above relation one arrives at 

JH /{pulse - cm) = 6.211 Joules. (6.21) 
The total power generated in the neck-down section of the horn alone is estimated to be: 

P neck = 2.5{Hz) ■ 6.211J/ (cm-pulse) ■ 86 cm = 1.335 kW. (6.22) 

The temperature rise in the horn conductor per pulse is estimated by calculating the 
volume per cm-length of the conductor, 

V = 2nr 5 = 2tt ■ 0.95 ■ 0.06525 = 0.38947 cm 3 . (6.23) 

The mass of aluminum in the neck-down volume (with density pal = 2.77 g/cc) is 1.0788 
g. The temperature rise in the horn conductor is estimated from, 

Q = mCpAT, (6.24) 

where, c p = 0.967 J/g-C. This leads to an estimated temperature rise per pulse AT = 5.95°C. 
As will be shown in the next section, the calculated temperature rise per pulse using a detailed 
finite element analysis is of the order of ~ 6°C. 

HORN Thermal Finite Element Analysis 



By focusing on the neck-down section, it is estimated that approximately 1.4 kW are 
generated due to joule heating. The temperature rise per pulse at the inner surface of the 
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conductor (based on aluminum density of 2.77 g/cc and c p = 0.967 J/g-C) is about ~6°C per 
pulse while the melting temperature of aluminum is 660 C. Numerical calculations based on a 
sophisticated finite element analysis verified the temperature rise in the inner horn conductor. 
Figures 16.151 - I6.19l depict temperature distributions and transients throughout the horn. The 
transient analysis of the integrated system is performed while considering all the participating 
mechanisms of heat transfer in order to find out the operating or "steady-state" temperatures 
of the system. As seen from the following results, the goal of maintaining the aluminum 
temperature below 100° C is achievable. It should be noted that with the replacement of 
RAW for the inner surface horn cooling with sub-cooled He the threshold of 100° C can be 
increased to several hundred degrees and thus reducing the demand for massive cooling since 
the heat transfer effectiveness will increase with higher allowable surface temperatures in the 
aluminum. 

Exploration of Alternative Horn Materials 



While aluminum is the baseline material for this study, attention has been focused to a 
new material called AlBeMet which is considered to be the highest performance aluminum 
alloy in the industry. The Figures 16.131 - 16. 141 provide its composition and properties relative 
to a number of other materials. In this on-going study this material will be furthered assessed 
for potential use in the horn structure. 

AlBeMet® 

AM162 

By weight, contains 62% commercially purs beryllium and 33% 
eommcrcraify pure aluminum 

Sy Metall u rgical ■deWnlti on, AlBeMet? is a n a Hoy but tan be 
considered £ composite 

i.'.r}- [;,<•■.'. s:-ieet: plate and bar are powder metallurgy products 

Th? pander i& produced tri a yy* vfvn'ilzatoii vrwwi wind' y-H'Iffo- yn^*ri^| 
p<?**)&r wiih fi fi-fly bfliyl!iiinn BtnJGLuto 

Thr pooler is d*i!5-1ifld enra* CjQnsalidatian pr-gca?±e£, «ach re*nlt3A$ in 
ditlip «ml mtClNnlcal pnptRifts. *hih) nurturing MBMItt* unlqg» pihyaieal 

Figure 6.13: Composition of AlBeMet alloy. 



Work in Progress 

In the course of this on-going study a number of issues still remain to be addressed in the 

effort to optimize the horn/target system and maximize its useful life. The "open" issues that 

are currently being addressed through material selection, detailed modeling and analysis and 

envisioned R&D efforts are summarized below. 
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AlBeMet® Property Comparison ] 
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Figure 6.14: Properties of AlBeMet Alloy. 




Figure 6.15: Geometry of the upstream end of the horn/target arrangement and temperature 
distribution in the narrow inner horn conductor prior to the arrival of the proton beam due 
to 240 kA current pulse. 
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Figure 6.16: Calculated temperature rise in the conductor skin depth (narrowest section) due 
to the half-sine 240 kA peak current pulse using finite element analysis. 



Transient temperatures over the two surfaces of the inner horn conductor 
T target = surface on the target side ; T_i = surface on current side 




Figure 6.17: Transient temperatures over the two surfaces of the narrowest section of the 
inner horn conductor. 
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Transient temperatures over the inner surface (current side) 
□f the horn at the front transition section 
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Figure 6.18: Transient temperatures over the inner surface of the upstream transition section 
of the horn. 



Transient temperature over the inner surface (current side) of the outer horn conductor 
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Figure 6.19: Transient temperatures estimated over the inner surface of the outer horn 
conductor. 
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• Gamma-Ray Heating of the Horn (especially of the inner conductor) Independent 
hadronic analysis are being performed to confirm the energy deposited in the inner 
conductor of the horn closest to the target. 

• Horn Thermal Stresses 

• Horn Magnetic Forces and Induced Stresses 

• Horn Vibration - Modal Analysis 

Detailed finite element analysis are under way to estimate the level of stresses generated 
in the horn, its modes of excitation and the means to optimize the horn structurally so it can 
meet all the physics criteria. 

• Horn Material (Aluminum) Irradiation Sc Fatigue It is envisioned that as part 
of the R&D and fatigue estimation effort, experience data will be reviewed further and, 
possibly, exposure of the material to dpa values will take place. Further, tapping on the 
expertise of collaboration members, the effects of the water environment on the changes 
induced on the aluminum surface and the issue of crack formation will be examined. 

• Discussion of Issues Related to Choice of Coolant In connection with the previous 
bullet, the effects and heat removal capacity of the choice of coolant fluids will be further 

assessed. 
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6.7 Conceptual Design of the Horn Power Supply 



The main objectives of the horn [20] power supply design is to achieve an output pulse am- 
plitude of 250KA, a pulse flat top of 2.5 /iS or longer with a 1% flatness and pulse to pulse 
repeatability, and a repetition rate of 2.5 pulses per second. 

Most commonly used scheme is the capacitor discharge type [21] [22] [22] [23] • In this type of 
circuit, a capacitor bank stores the energy, and a main discharge switch release the energy to 
the load through transmission lines. For very long distant transmission, pulsed transformers 
have been added into the KEK design and CNGS horn system. The horn as an electrical 
load is usually being described as an inductor in series with a resistor. Discrete parameters 
of inductance and resistance are also used to formulate the short length, low impedance 
transmission lines when associated with low bandwidth pulse. Hence, the circuit can be 
simplified clS db RLC discharge circuit. 

The pulse rise time, T r , is usually approximated by the quarter period of sine wave. For 
a lossless LC oscillation circuit, this can be determined by the equation T r = | \J LC , and 
the load current is given by 

/ W = ^\/ / f sin(-=t). (6.25) 



The maximum current amplitude is 7 max = V -y/^for lossless LC circuit. 




In the case of critical damped RLC circuit, R = 2 ^JL/C and 



I(t) = e~ aK (6.26) 




The maximum output current can be determined by 

= Wi e_1 ' (6,2T) 

and the current rise time is T r = ^ = 

Compare two cases, for a given L, C and charging voltage V D . the output current of the 
critical damped one has peak amplitude of 0.3679 of the lossless one. To reach the same 
maximum output current, the initial voltage of the critical damped circuit has to be 2.7183 
times of the lossless one. It implies that, for the critical damped one, the energy stored in the 
capacitor bank as well as the charging power has to be 7.389 times larger. 

In the high current path, the resistance caused voltage drop and energy dissipation are 
critical factors to be considered. The cooling system for heat removal from the effective 
resistor, and the additional charging power required to make up the resistive dissipation can 
be very costly. Hence the low circuit resistance is preferred. 



6.7.1 Effective Resistance & Skin Effect 

The load and transmission line resistance varies with frequency due to skin effect. For any 
given material, the skin depth S s is inverse proportional to the square root of the frequency f, 
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and the effective resistance R e // is proportional to the oscillating frequency. Here 

S s = , , l (6.28) 

and 

R cff = (6.29) 

where 1 is the conductor length and b is the conductor width. The material's conductiv- 
ity, relative permeability, and the free space permeability are described by a, jiR, and /i , 
respectively. 

One can see that slower frequency leads to lower effective resistance. The other factors 
associated with effective resistance are the length, width, and permeability and material con- 
ductivity. For non-magnetic material, the relative permeability is close to unit. The switching 
device on-state resistance and hardware connection joints resistance also contribute to the 
total resistance. 

In summary, the lower effective resistance can be achieved by using: 

1. Lower frequency, higher skin depth; 

2. Wider conductor width; 

3. Shorter conductor length; 

4. Higher conductivity material; 

5. Lower switch on-state resistance; and 

6. Lower connection joint resistance. 

6.7.2 Inductance Issues 

The total inductance includes the horn inductance, transmission line inductance, series in- 
ductance of the capacitor, switch inductance, and circuit loop stray inductance. The external 
inductance depends on the inductor geometry and material permeability. The internal induc- 

dl(t) 



dt 



tance has frequency dependence. The voltage across the inductor is defined as Vz(t) = L 
Hence, the larger inductor and faster current rate of change requires higher voltage. 

The current going through the inductor is the inductively stored energy. For a lossless 
circuit, the total inductive energy shall be equal to the total capacitive energy storage, i.e. 

\LI 2 = \CV 2 . (6.30) 

Therefore, we have that the higher the inductance and current, the higher the capacitance 
and its initial voltage. 

For a reasonable design, the total inductance shall be kept as low as possible, and the 

current rise time shall be chosen to accommodate the device operating voltage. 
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6.7.3 Principle Design Example 

Let us consider a system with overall inductance of 2.5 /iH, a total resistance of 2 mf2, and a 
capacitor bank of 16 mF. If the initial capacitor voltage is 4000 Volts, the peak output current 
amplitude is about 288 kA with 314 /j,s rise time. 



V 



m 




L1 



Figure 6.20: Simplified circuit diagram. 



The stored energy in the capacitor is 128 kJ. For 2.5 Hz pulse repetition rate, the minimum 
charging power supply is 320 kW. 

For the 3000-series aluminum being considered in the Horn mechanical design, the material 
skin depth is 0.06525 cm at 25 kHz. Derive from it, we have the skin depth is 0.3574 cm at 
833 Hz. 

The parameters used in this example are rough estimates based on similar systems, as 
listed on Table 16.31 The resistance used in the example is tight for the chosen frequency. If 
higher resistance has to be used, then the voltage and capacitance have to be increased. 



6.7.4 Major Components 

The major components of this system include: 

1. Charging Power Supply; 

2. Capacitor Bank; 



3. Discharge Switch; 
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Figure 6.21: Circuit simulation with V=4000, Cl=16 mF, Rl= 2 m SI, Ll=2.5 /iH, L2=0.1 
/iH, R2=7.5 ma 



4. Reverse Diode; 

5. Transmission Line. 

In this design, Self-Healing capacitors are being considered for the fault tolerance, and 
increased reliability. This type of capacitor is usually rated under a few kilo-volts. Multiple 
capacitor cells have to be used to divide capacitor bank into smaller units with lower stored 
energy per cell for safety concerns. 

The discharge switch in favor is the light triggered thyristor. The newest in this category 
is the EUPEC T2563N80. The advantage of this device is its high voltage rating of 8000 V, 
high forward current rating of 5600 A rms, and 63 kA surge current. Multiple thyristors have 
to be used in parallel to carry the 250 kA pulse current. 

This thyristor features light triggered gate structure. It eliminates the high voltage iso- 
lation trigger transformers normally used in conventional thyristors, and improves the high 
voltage hold-off and the noise immunity. Large reverse diodes are available from the same 
company. 

The trend of new designs is to use solid-state switch, which has much longer lifetime com- 
pared to gas discharge switches. The very high current capability required in this system limits 
the selection to thyristor types for its high power rating per single unit and cost effectiveness. 

The traditionally used ignitron is mercury vapor filled device. With rapid advancement of 
solid state devices, it is being replaced by thyristors. Other solid state devices, such as IGBT 
and MOSFET are limited by their current capability. 

The transmission line being considered is similar to FERMI NuMI j2Hj design for its ultra 
low resistance. 
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The design options of high voltage, high current pulsed system are often limited by the 
industry development and available components. In this case, the preferred operating voltage 
is under 5 kV. The total resistance shall be kept to less than or around 2 mQ. 
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7 Conventional Facilities and Target Hill 

The proton beam will be extracted from the AGS and will utilize part of the RHIC beam 
transport before exiting the decommissioned neutrino beam line tunnel in a northerly direction 
and at an upward angle of approximately 13.8 degrees. The beam will bend towards the west 
by approximately 68.5 degrees, then down a total of 25.1 degrees to the proton target. A 200 
m decay region follows with the neutrino beam entering the hadron beam stop at 11.26 degrees 
to the horizontal. The Near Detector Facility is located 285 m from the target, 21 m below 
ground level. A plan view of this arrangement at BNL is illustrated in the Figure I5~2k nd the 
overall beam facility is shown in Figure 17.11 

The vertical beam geometry results in a hill 50 meters high at the apex of the proton 
beam. This geometry provides for the neutrino beam's 11.3-degree entry into the earth and 
avoids potential irradiation of soil close to the Long Island water table. This, together with 
impermeable rain water barriers to prevent rain water penetration of potentially irradiated 
soil, is consistent with present ground water protection practices at BNL. A simplified vertical 
beam profile and a cross-section of the hill is shown in Figures 15.11 and 15.31 




Power supply • ••" Pion Decay Pipe 

Building . /'\ 



Beam Stop 

^S^ppa^ ^ near Detector 

< T NORTH 



Figure 7.1: The hill, truncated ~2 meters below the beam and the new beam transport 
elements are illustrated. 



The new conventional facilities can be subdivided into the following areas: 

• Proton beam transport 

• Target area 

• Decay tunnel 
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Figure 7.2: 13.8 Degrees vertical bend. 




Figure 7.3: 68.2 Degrees horizontal bend. 




Figure 7.4: 16.6 Degrees vertical bend. 
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• Hadron beam stop 

• Near Detector Facility 



7.1 Proton Beam Transport 

Existing utilities and roads will be relocated and approximately 726,350 cubic meters of sand 
fill will be placed forming the hill. The sand fill will be placed in 0.3 m lifts and compacted to 
98% of its maximum density. The fill will be placed early in the project allowing it to settle 
for several years before re-excavation for placement of the tunnel. Approximately 330 m of 
3-m diameter tunnel, overburdened with 6.0 m of fill will be required for proton transport. 
There will be a waterproof liner installed .6 m below the surface of the overburden. Portions 
of beam transport are depicted in Figures I7.2| 17.3^ and 17.41 

Access to the tunnel will be provided at the beginning of the vertical rise and near the 
target area. Equipment transport up & down the slope will be facilitated by use of a motorized 
railway on the aisle side of the transport elements in the tunnel. A typical section of the proton 
tunnel is shown Figure 17.51 
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Figure 7.5: Typical proton beam tunnel section. 



Two power supply/utility buildings will be provided, one located low near the existing 
beam line, the other located high near the target area. These buildings will house power 
distribution systems, power supplies, water pumping systems, instrumentation and controls 
for the beam line. 

The cooling water system will use a 3.5 MW cooling tower for primary heat rejection with 
four isolated, closed loop cooling systems for: 

• All transport magnets (~2.8 MW) 

• Two power supply areas (~. 5 MW) 

• Horn cooling (~. 2 MW) 
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Each system will contain redundant pumps, a heat exchanger, a full flow filter and a 
side stream deionizer capable of maintaining the system at 2-5 Megohm-cm. The system 
controls will be PLC based and be capable of monitoring and reacting to water leaks if 
they occur. In general, electrical power will be distributed around the site at 13.8 kV. Unit 
substations will transform the power into convenient voltages, typically 480 and 208/120 
volts. Electrical power is divided into two major categories: conventional and experimental. 
Conventional power encompasses building power for lighting and convenience power, HVAC, 
and miscellaneous equipment; and site distribution. Emergency power will be provided as 
required from existing circuits. 

Experimental power feeds all the power supplies and associated equipment, such as, cool- 
ing water pumps, cooling towers, etc. Designs will follow the requirements of the National 
Electrical Code and industry standards. 



Table 7.1: Preliminary power requirements. 



Location 


Conventional 


Experimental 


Total 




Power (KVA) 


Power (KVA) 


Power 
(KVA) 


Lower Power House 


150 


2750 


2900 


Upper Power House 


150 


2050 


2200 


Near Detector 


100 


900 


1000 


Total Power (KVA) 


400 


5700 


6100 



The present radiation security system in the RHIC transport line will be extended and 
three new radiation security access gates will be installed. Access to the new tunnel will be 
allowed during RHIC operations by implementing the necessary features such as shield walls 
and beam plugs in the existing tunnel. A fence will be provided to limit access over the beam 
tunnel while the beam is in operation. 

Fire detection and protection will be provided where appropriate throughout the facility 
and connected to the BNL site fire alarm system. High sensitivity smoke detection systems 
will be utilized in increased risk or high value 

The proton beam tunnel will be provided with dehumidification, ventilation and a com- 
pressed air system for device controllers. Heat and air-conditioning will be provided in the 
power supply buildings. All buildings will be connected to the BNL site sanitary system where 
the appropriate monitoring of effluent is provided. 

Ground water monitoring wells will be provided to insure compliance with all Local, State 
and Federal ground water protection requirements. 

7.2 Target Area 

The target area will include the proton target, horns, horn power supply and water cooling 
system. The facility will be a high bay building with a stepped foundation and an overhead 
crane of sufficient capacity to handle the largest load of approximately 40 tons, see Figure 
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17.61 A shielded storage area will be provided for radioactive component storage and repair. 
Modular concrete and steel shielding will provide radiation shielding. Access to the horn 
vault for installation and removal of the horns is accomplished by removing the modular 
shielding. Present plans call for hanging the horns from shielded supports with the ability to 
survey and connect or disconnect the horns from above the shield in a relatively low residual 
radiation environment as shown in Figure 17.71 The design of this SX6£L £lS well as all 
will incorporate ALARA radiation principles. A collimator downstream of Horn j^l will be 
installed to intercept a portion of the off-axis beam that would otherwise interact in the soil 
along the decay tunnel. 



5. v 




TARGET AREA 

Figure 7.6: Target area. 

Survey tolerances, both absolute and relative, are driven by the physics requirements which 
are less stringent than for the "MINOS" experiment in Soudan, MN.. The geodetic position 
of BNL is well known and tied to the National Geodetic Survey system. The geodetic position 
of the beam front-end elements and the Detector site will be taken in 2 to 3 simultaneous 
epochs giving the relative positions of the two sites within a few centimeters. A more difficult 
issue will be in determining the position of the subterranean detector relative to the surface. 
Gyro-theodolite observations coupled with traditional line-of-site observations (if possible) will 
provide the data necessary. If line-of-site observations are not possible, an Inertial Navigation 
System may be used giving the relative position within a few meters. 
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Figure 7.7: Cross section at horn#2. 



7.3 Decay Tunnel 

The decay tunnel will be a 2-meter diameter steel tube; 185 meters long, with seal welded 
joints. There will be a thin helium window at the downstream end of Horn #2 and a helium 
window at the upstream end of the hadron beam stop. Simulations will be done to predict the 
temperatures and stresses on the windows and the appropriate measures taken to insure they 
withstand the radioactive environment. The contained volume will be purged with helium 
gas. Access to the upstream window will be provided through the target area. There will 
be no utilities or access to the decay tunnel between helium windows. The tunnel will be 
overburdened with 9 meters of earth fill with a waterproof liner installed approximately 0.6 
meters below the surface, extending 9 meters on each side of the tunnel. 

7.4 Hadron Beam Stop 

The hadron beam stop will be approximately 7 m wide x 7 m high x 9 m long. Elevation 
view is shown if Figure 17.31 It will have a stepped foundation to approximate the 11.3- 
degree downward angle. The lowest portion of the beam stop will be approximately 3m 
above the ground water table. The stop will consist mainly of existing steel plate from the 
decommissioned neutrino beam line at BNL, overburdened with 4m of soil. A waterproof 
liner will be placed .6m below the surface of the soil to insure no rainwater can penetrate. 
Simulations will be carried out to predict the temperatures and stresses on the front face of 
the stop and the appropriate materials chosen. It is presently assumed that no cooling of the 
material is necessary and only temperature interlocks are required. 

7.5 Near Detector Facility 

The near detector facility will be located 285 m from the target, 21 m below ground level. 
The near detector facility is depicted in elevation view in Figure 17.91 The facility will consist 
of welded steel tunnel sections 6 m in diameter, with a 5 m access shaft to the surface. A 9 m 
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BEAM STOP AND DECAY TUNNEL 



Figure 7.8: Beam Stop and decay tunnel. 

wide x 15 m long service building will be constructed over the access shaft with a removable 
roof to access equipment below with a mobile crane. This building will contain electrical 
distribution, HVAC units, water cooling systems, power supplies and experimental equipment 
associated with the detector. An elevator will be provided for accessing the detector from the 
service building. A unit substation and cooling tower will be provided for utilities in this area. 
Estimated requirements are: 

Power 1000 KVA 

Cooling - - 500 kW 

Since the facility is located below the water table, installation of the tunnel sections will 
involve excavation of soil to the water table, installation of sheet piling, excavation and de- 
watering of the sheet pile site, and installation of the tunnel sections. The tunnel is then seal 
welded and back-filled with soil. 

7.6 Radiation Shielding for the Super Neutrino Beam 

The conceptual shielding design for the Super Neutrino Beam meets the necessary standards 
for radiation protection and provides a basis for the initial cost estimates for the project. The 
design has been developed using existing reports, calculations, and simple analytic techniques 
to scale to the proposed facilities for the proton transport, target station, decay tunnel, and 
beam stop. More detailed modeling will be conducted when the design of the various compo- 
nents are more mature. The shielding has been designed to: 

• Meet standards for chronic exposure to people at adjacent areas and off-site. 



Prevent contamination of the ground water. 
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NEAR DETECTOR 
VERT ICAL SECTION 



Figure 7.9: Near detector vertical section. 

• Reduce direct exposure sufficiently so that personnel can service equipment at the facility 
during beam operations. 

A brief description of the methods and assumptions used will be given in this section. 
7.6.1 Source Terms 

The shielding design requires estimations of the source terms and their locations. Logically, 

the Super Neutrino Beam breaks into four simple areas, proton transport from the AGS to the 

target station, target station, decay tunnel, and the beam stop. The proton beam transport 

is expected to have an integrated loss of 0.1% with less than 0.03% lost in any one location. 

The target station has been designed assuming 100% of the beam interacts in the target. It is 

assumed that no beam interacts with the helium in the decay tunnel. In reality, approximately 

1% of the incident beam will survive the target and interact in the helium. These neglected 

interactions are compensated for by the assumption that 100% of the beam interacts in the 

target. 20% of the beam will survive the target and strike the beam dump. 

An MCNPX calculation by A. Stevens has been used to estimate the initial radiation 

pattern in the near shield. The calculation assumed that the 25 GeV proton beam struck an 

iron target 1-meter long and 3 cm in radius. This target was placed in a tunnel with a 1.5 

meter radius and surrounded by soil. The initial radiation pattern has then been extrapolated 

using the analytic techniques of Teschj2Z] and discussed by A.H. Sullivan|2*S]. This method 

was used for the beam transport, target station, and decay tunnel. 
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The beam stop was designed by using an existing design [21] of a 29 GeV beam dump in the 
AGS slow beam area. The analytic techniques were then used to extrapolated were necessary. 

7.6.2 Ground Water Protection 

The interaction of high-energy protons will create radioactive products in the surrounding 
shield. It is important that these radioactive products are not transported to the water 
table. The BNL standardjnOl requires that the design of this facility prevent radioactive 
products in the groundwater from exceeding 5% of the drinking water standard. The use 
of modular concrete and steel shielding for the entire project is cost prohibitive. Most of 
the facility will have an earthen shield. Modular shielding will be used at the locations of 
highest loss and when possible to reduce the amount of radioactive products created in the 
surrounding soil. The design standard allows for the use of approved capping materials to 
prevent created radioactive products from being transported by rainwater to the water table. 
Caps constructed of geomembranes or concrete will be used to meet the BNL standard for 
groundwater protection. 

22 Na is the radioactive isotope that will determine the design of the cap which protects the 
groundwater. Based on the transport model used in the BNL standard, the cap must prevent 
rainwater from leaching the soil where the created 22 Na in a year exceeds 1.05*10 6 22 Na 
atoms/cc. This corresponds to a neutron flux (E>20 MeV) of 9.4*10 -13 neutrons/cm 2 per 
incident proton. 

The loss assumption for the proton transport has a maximum local chronic loss of 0.03% 
of the proton beam. The cap must have at least 4.6 meters of soil between it and the tunnel 
wall. For the present conceptual design the membrane has been placed at 5.5 meters from the 
tunnel wall giving a factor of 5 more reduction then required by the standard. The membrane 
is covered with 0.6 meters of earth. 

It is assumed that 100% of the beam interacts at the target station. If it were sited inside 
a tunnel shielded by dirt the liner would need to be greater than 9 meters from the target. 
The total amount of activation in the soil would be substantial. It was decided that the target 
area should have the shield constructed from modular steel and concrete blocks. The target 
area shield will have a thickness equivalent to 0.6 m of steel and 4.3 m of heavy concrete. The 
target area will be inside a building. The concrete floor will prevent rainwater from leaching 
activation products from the soil underneath the target station. The roof of the target building 
will prevent rainwater from entering the concrete and steel shielding. Additional shielding will 
be placed underneath the target and horns to reduce the activation of soil underneath the 
concrete floor as shown in Figure 17.61 

The decay tunnel has been designed with the membrane 8.5 meters from the tunnel wall. 
The membrane will have 0.6 meters of soil overburden. At this stage of the design it was 
decide to be conservative in the estimate for the decay region. It is expected when the final 
design of the target area is complete that the membrane can be closer to the tunnel and the 
berm thickness decreased. 

The beam dump is constructed from plates of steel. The dimensions of the steel were 
chosen to reduce the neutron flux sufficiently such that the activation of the soil in contact 
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with the edges of the steel satisfies the BNL standard. This requires the iron beam dump to 
have a radius of 3.4 meters and a length of 9.2 meters. The water cap for the decay tunnel 
will be extended over the beam dump, as shown in Figure 17.31 This is necessary to prevent 
water from percolating through the beam dump. 

There are four monitoring wells planned for the Super Neutrino Beam area. The beam 
dump, decay tunnel, target station, and proton transport will each have an appropriately sited 
monitoring well. In addition, existing wells will monitor the portion of the proton transport, 
which is close to the AGS. 

7.6.3 Sky Shine 

Dose to adjacent facilities occurs due to leakage of neutrons through the shielding. The dose 
to people at adjacent facilities is calculated using the sky shine formulation discussed by A.H. 
Sullivan [2*%]. The exposure at a location more than 100 meters from the source is given by: 

H = 70(H A) exp(-i?/600)/(i? 2 ), (7.1) 

where H is the dose rate in mrem/hr, Ho A is the integrated dose through the berm in rem-m 2 
per hour, and R is the distance from the source to the occupied location in meters. 

The limits are established by integrating the exposure over a year of operation. In addition 
occupancy factors can be used were appropriate. Five locations were examined for the yearly 
exposure from the super neutrino beam. The assumptions on occupancy and exposure limits 
are given in the Table 17.21 



Table 7.2: Occupancy and exposure limits. 



Location 


Occupancy 


Occupancy per 


Facility limit 




hours/ week 


work hour 


mrem / yr 


Off-Site 


168 


1 


5 


Linac 


40 


1 


100 


Blip 


40 


0.25 


25 


BAF 


40 


0.25 


25 


Bldg. 1005 


40 


1 


25 



Based on these occupancy assumptions the exposure limit for off-site is the most stringent 
for all loss locations. The target station dominates the contribution to off-site dose and is less 
than 10% of the allowed limit. 

7.6.4 Direct Exposure 

Dose rates external to the shielding will typically be less than a few mrem/hr under normal 
operating conditions. The highest localized dose rates occur at the target building and are 
less than 60 mrem/hr . Barriers, postings and procedures will control access to the shielding. 
Abnormal operating conditions can cause substantial increases in the dose rates. Elements, 
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which can cause abnormal conditions, will be monitored and a fast beam interrupt will limit 
the duration of the abnormal condition. Radiation detectors will be distributed along the 
facility to prevent personnel exposure and limit fault conditions. A dual PLC based access 
control system will prevent access to the transport tunnel, except under allowed conditions. 
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8 Cost Estimate and Schedule 

The cost estimate of the AGS-Based super Neutrino beam Facility has been performed fol- 
lowing the "bottoms up" approach. After the performance goals and physics design of each 
technical systems are completed, cognizant engineers who have built similar systems were as- 
signed to prepare a cost estimate of sufficient detail and based on previous experience. Each 
estimate typically includes the cost of the detailed engineering design, procurement, manu- 
facturing, testing and installation. Material and labor costs are captured as separate entries 
for each process. 

Most of the cost numbers are based on BNL's recent experience in building RHIC, SNS 
ring, and LHC magnets. They also draw on extensive searches of price given in catalogs 
and on vendor's quotations. Although, a formal Work Breakdown Structure (WBS) has not 
been implemented yet, all systems have been estimated by going down to the component and 
assembly level. From our past experience, both the manufacturing approaches and the cost 
estimates are consistent with good engineering practices and are credible. 

The costs given reflect the 'direct cost" of the item. In other words, estimates do not 
include G & A (lab overhead), contingency and ED &I. These cost elements are added for the 
direct cost following guidelines using standard estimating. 

All the figures are in FY 2003 US Dollars. No inflation is included since the construction 
period is not known at this time. 

A preliminary estimate of the direct cost of the accelerator system is given Table 18.11 



Table 8.1: Accelerator system cost. 



1.2 GeV SC Linac: 




Front End 


$2.5 M 


LE SC Linac 


$38.3 M 


ME SC Linac 


$30.7 M 


HE SC Linac 


$28.1 M 


AGS Upgrades: 




AGS Power Supply 


$32.0 M 


AGS RF Upgrade 


$8.6 M 


AGS Injection Channel 


$3.7 M 


Full Turn Extraction 


$5.5 M 


Shielding 


$3.2 M 


Installation 


$4.2 M 


Total Direct Cost 


$156.8 M 



The direct cost of the neutrino beam systems is summarized in Table 18.21 
The resultant total direct cost of the 1 MW AGS super Neutrino Beam facility, not in- 
cluding both near and main detectors, is about $ 218.5 M. More detailed cost estimate for 
each system are given in Appendix B. The preliminary total estimated cost (TEC) is $369 
M in FY03 dollars, including EDIA 15 %; contingency 30 %; BNL project overhead 13 %. 

Escalation cannot be estimated without a project start year. 
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Table 8.2: Estimated neutrino beam cost. 



Item 


Basis 


Cost 


rroton transport 


T ) T T T /" i J • i 

RHIC Injector 


$14.8 M 


target /Horn 


E889 


(T 1 r r TV /T 

$5.5 M 


Shielding/Dump 


1\ cw 


(T? r o A vr 

$5.o M 


Decay Tunnel 


E889 


$0.4 M 


Hill. Const. 


New 


$13.0 M 


Near Detector Vault 


RHIC 


$7.3 M 


Conventional Facil. 


RHIC 


$7.5 M 


Other Const. 


New 


$2.2 M 


Installation 


RHIC 


$5.2 M 


Total Direct Cost 




$61.7 M 



It is estimated that three years of R & D are needed to build prototypes and perform detail 
engineering designs to reduce cost and improve operation reliability. This will be followed by 
4.5 years of construction and 0.5 year of commissioning to get this facility ready for physics 
research. The planned schedule is represented in Table 18.31 Where FY1 represents the first 
year of project approval. 



Table 8.3: Schedule. 





R&D 


Construction 


Research 




FY1 


2 


3 


4 


5 


6 


7 


8 


9- 15 


R&D 


XX 


XX 


XX 














Construction 








XX 


XX 


XX 


XX 


X 




Commissioning 
















X 




Research 


















xxxx 
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A Appendix A: Design Parameters 
A.l Facility Level Parameters 



Table A. 4: Facility level parameters. 



Proton Beam Energy 


28 GeV 


Protons per Pulse 


8.9xl0 13 


Average Beam Current 


37.5 uA 

r 


Repetition Rate 


2.5 Hz 


Pulse Length 


2.58 usee 

r 


Number of Bunches 


23 


Protons per Bunch 


3.87xl0 12 


AGS Circumference 


807.1 m 


Bunch Length 


40 nsec 


Bunch Spacing 


60 nsec 


Extraction Gap 


160 nsec 


Average Beam Power 


1 MW 


Normalized Emittance-X 


100 7r mm-mrad 


Normalized Emittance-Y 


100 7r mm-mrad 


Longitudinal Emittance 


5.0 eV-sec 


Energy Spread (AE/E) 


0.001 


Target Material 


carbon-carbon 


Target Radius 


6.0 mm 


Target Length 


80 cm 


Beam Size on Target 


2 mm (rms) 


Beam Elevation at Target 


43 m 


Decay Tunnel Length 


200 m 


Beam Dump Length 


20 m 


Distance of Near Detector from Target 


275 m 


Physics Time per Year 


10 7 sec /year 


Number of Protons on Target 


2.2 xlO 21 p/year 
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A. 2 Front End and Warm Linac Parameters 



Table A. 5: Front end and warm linac parameters. 



Energy of Warm Linac 


200 MeV 


Normalized Horizontal Emittance (rms) 


2.0 7r mm-mrad 


Normalized Vertical Emittance (rms) 


2.0 7r mm-mrad 


Longitudinal Emittance (rms) 


0.125 MeV-deg 


Macro-Pulse Average Current 


21 mA 


Macro-Pulse Peak Current 


28 mA 


Repetition Rate 


2.5 Hz 


Pulse Length 


720 usee 


Chopping Rate 


0.75 % 


Energy Spread (AE/E) 


0.002 


Energy Jitter (SE/E) 


0.001 
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A. 3 Superconducting Linac Parameters 



Table A. 6: Superconducting linac parameters. 



Linac Section 


LE 


ME 


HE 


Kinetic Energy Initial/Final, MeV 


200/400 


400/800 


800/1200 


Frequency, MHz 


805 


1610 


1610 


No. of Protons / Bunch x 10 8 


8.70 


8.70 


8.70 


Temperature, °K 


2.1 


2.1 


2.1 


Cells / Cavity 


8 


8 


8 


Cavities / Cryo- Module 


4 


4 


4 


Cavity Separation, cm 


32.0 


16.0 


16.0 


Cold- Warm Transition, cm 


30 


30 


30 


Cavity Internal Diameter, cm 


10 


5 


5 


Cell Reference /3q 


0.615 


0.755 


0.851 


Cell Length, cm 


11.45 


7.03 


7.92 


Length of Warm Insertion, m 


1.079 


1.379 


1.379 


Total No. of Periods 


6 


9 


8 


Length of a period, m 


6.304 


4.708 


4.994 


Total Length, m 


37.82 


42.38 


39.96 


Cavities / Klystron 


1 


1 


1 


No. of rf Couplers / Cavity 


1 


1 


1 


Coupler rf Power, kW (*) 


263 


351 


395 


Total No. of Klystrons 


24 


36 


32 


ZoTq, ohm/m 


378.2 


570.0 


724.2 


Qo x 10 10 


0.97 


0.57 


0.64 


Ave. Axial Field, E a , MV/m 


13.4 


29.1 


29.0 


Ave. Dissipated Power, W 


2 


11 


8 


Ave. HOM-Power, W 


0.2 


0.5 


0.4 


Ave. Cryogenic Power, W 


65 


42 


38 


Total Ave. RF Power, kW (*) 


17 


31 


30 


Norm, rms Emittance, ir mm-mrad 


2.0 


2.0 


2.0 


Rms Bunch Area, vr °MeV (805 MHz) 


0.5 


0.5 


0.5 



(*) Including 50% RF power contingency. 
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A. 4 AGS Parameters 



Table A. 7: AGS parameters. 



Injection Energy 


1.2 GeV 


Extraction Energy 


28 GeV 


Number of Injection Turns 


240 


Stripping Efficiency 


98% 


Beam Bump 


25 kW 


Electron Collection 


25 Watts 


Number of Protons 


8.9 xlO 13 


Harmonic Number 


24 


Filled Bucket 


23 


Repetition Rate 


2.5 Hz 


Strip Foil 


carbon-carbon 300 /xg /cm 2 


RF Frequency at Injection 


8.0 MHz 


RF Frequency at Extraction 


12.0 MHz 


Peak RF Voltage 


1.0 MV 


Transition Gamma 


8.5 


Horizontal Tune 


8.7 


Vertical Tune 


8.9 


Beta Maximum 


22.5 m 


Beta Minimum 


11.5 m 


Horizontal Aperture 


7 cm 


Vertical Aperture 


7 cm 


Power Supply 


4300 A, 25 kV 
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A. 5 Target/Horn Parameters 



Table A. 8: Target/Horn parameters. 



Target Material 


carbon-carbon 


Target Radius 


6.0 mm 


Tarsret Lenfrth 


80 cm 


Heat Deposition 


7 3 k.T / milse 


Coolinff Medium 


He Gas 


Operating Peak Temperature 


800°C 


Horn Small Radius 


70 mm 


Horn Large Radius 


610 mm 


Horn Thickness 


2.5 mm 


Horn Length 


217 cm 


Horn Current 


250 kA 


Horn Peak Current Density 


690 kA/cm 2 


Repetition Rate 


2.5 Hz 


Power Supply Wave Form 


Sinusoidal, base width 1.20 msec 


Power Supply Voltage 


5 kV 


Average Power 


700 kW 


Heat Deposition on Horn 


3 kJ/Pulse 


Cooling Medium of Horn 


Water 
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A. 6 Decay Tunnel and Shielding Parameters 



Table A. 9: The Decay tunnel and shielding parameters. 



Decay Tunnel Radius 


1 m 


Decay Tunnel Length 


200 m 


Beam Dump Material 


Steel 


Beam Dump Length 


9 m 


Shielding of Target Bldg. 


0.6 m Steel and 4.3 m of 
Heavy Concrete 


Shielding of Decay Tunnel 


9 m Soil 
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A. 7 Conventional Facilities and Target Hill Parameters 



Table A. 10: Conventional facilities and target hill parameters. 



The Linac Tunnel Length 


120 m 


Linac Tunnel Diameter 


3 m 


Klystron Gallery Length 


120 m 


Klystron Gallery Dimension 


22mxl2m(W*H) 


SRF Testing and Assembly Building 


12 m x 36 m 


Shielding of Linac Tunnel 


3 m Soil 


Linac Service Building 


12 m x 22 m 


Beam Transport Tunnel 


330 m 


Beam Transport Tunnel Diameter 


3 m 


Shielding of Beam Transport Tunnel 


6 m Soil 


Beam Transport Service Buildings 


2x(250 m 2 ) 


Beam Elevation at Target 


43 m 


Target Hill Height 


52 m 


Target Hill Base Width 


154 m 


UP Hill Length to Target 


200 m 


Down Hill Length from Target 


220 m 


Target Service Building 


16.5 m x 15 m 


Electrical Substation 


6100 KVA 


Water Cooling Capacity 


3.5 MW 


Service Road 


250 m 
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B Appendix B: System Cost Estimate 
B.l Linac Upgrade Cost 



Table B.l: Linac upgrade cost. 



Direct Costs 


k$ 


Ionsource Power Supplies 


80 


Driver Stage RF Systems 


123 


RF Modulator PS's 


23 


High Power Cap Banks 


730 


Crowbar System 


76 


Low Level RF Amps 


102 


Quadrupole Pulsers 


1016 


Misc. Cabling, etc. 


27 


Labor 


325 


TOTAL DIRECT 


2,542 
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B.2 SCL Cost 



Table B.2: Cost ('02 $) and other parameters. 





LE-Section 


ME- & HE-Sections 


Frequency, MHz 


805 


1,610 


AC-to-RF Efficiency 


0.45 for Pulsed Mode 


Cryogenic Efficiency 


0.0014 @ 2.1 °K 


'FTIpptripitv (!n^t 


0.05 $/kWh 


Linac Availability 


75 % of yearly time 


Nnrinfl 1 (!nn nnpt (Irmt 

liUl 111CX1 V Willi \JL\j Li . vJUO Li 


20 k$/m 50 k$/quad 


^n'HPT'pnTirinpfi'nrf ( /~\G"t" 


200 k$/m 


100 k$/m 




360 k$/module 


200 k$/module 


C ,ct^it of Pvlvcrt"TOn in kS* 
vjuoi ui i vi^y o li vjn^ 111 i\.yJ . 






Pnw/pr Snnrpp 

1 UVV CI kJ \J LL1 l^C 


170 


160 


Tt ^ n s -1- PS 

XI ClillQ. | 1 kJ 


140 


130 


Clircul -r-Loads 


95 


80 


Controls 


70 


60 


Waveguides 


25 


20 


Total 


500/Unit 


450/Unit 


Cost of Refrigeration: 






Refrigeration 


8.0 M$ 




Control 


1.0 M$ 




Ancillary 


1.0 M$ 




Labor 


6.0 k$/m 




Transfer Lines 


18.0 k$/m 




Electrical Distribution 


0.14 $/W (AC) and 10 k$/m 


Control 


10% of the Total Cost above 


Tunnel Cost 


15 k$/m (Linac) and 40 k$/m (Klystr.) 
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Table B.3: Cost estimate of the SCL project. 



Linac Section 


LE 


ME 


HE 


Capital Cost '02 M$: 








RF Klystron 


12.000 


16.200 


14.400 


Electric Distr. (**) 


0.435 


0.438 


0.413 


Refrigeration Plant 


(+)10.908 


3.017 


2.959 


Warm Structure (**) 


0.733 


1.156 


1.040 


Cold Structure 


8.429 


4.796 


4.492 


Control System 


3.250 


2.561 


2.330 


Tunnel (**) 


2.560 


2.564 


2.417 


Total Cost, '02 M$ 


38.315 


30.731 


28.052 


Operation Cost, '02 k$/y (***) 


54.6 


64.9 


61.4 



(**) Includes 4.5m long matching insertion between LE and ME sections 
(***) Includes only electric bill for RF and Cryogenic Power converted to AC 
(+) Includes 8.0M$ for Main Cryogenic Station 



Table B.4: Level 3 of WBS and cost (in 2002 k$) of total SCL. 



The 1.2-GeV Linac 


Material 


Labor 


Total 


2.1.3 


Low-Energy SCL 


29,116 


9,199 


38,316 


2.1.4 


Medium-Energy SCL 


23,047 


7,685 


30,732 


2.1.5 


High-Energy SCL 


20,991 


7,060 


28,051 




Total 


73,154 


23,944 


97,099 
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Table B.5: Levels 4 and 5 of WBS and cost 



(in 2002 k$) of LE-section of SCL. 





Low-Energy SCL 




Material 


Labor 


Total 


2.1.3.1 


Cryo-Modules 




6,743 


1,686 


8,429 




2.1.3.1.1 


Tanks 


2,248 


562 


2,810 




2.1.3.1.2 


RF Cavities 


4,495 


1,124 


5,619 


2.1.3.2 


Refrigeration (*) 




8,726 


2,182 


10,908 


2.1.3.3 


Warm Insertions 




586 


147 


733 




2.1.3.2.1 


Quadrupoles 


352 


88 


440 




2.1.3.2.2 


Steering 


59 


15 


73 




2.1.3.2.3 


BPM 


59 


15 


73 




2.1.3.2.4 


Vacuum Sys. 


117 


29 


147 


2.1.3.4 


RF Power Supply 




9,600 


2,400 


12,000 




2.1.3.3.1 


Klystrons 


2,880 


720 


3,600 




2.1.3.3.2 


Transmitters 


2,304 


576 


2,880 




2.1.3.3.3 


Circulators 


1,824 


456 


2,280 




2.1.3.3.4 


Controls 


1,344 


336 


1,680 




2.1.3.3.5 


Waveguides 


480 


120 


600 




2.1.3.3.6 


Couplers 


768 


192 


960 


2.1.3.6 


Electrical Distribut. 




348 


87 


435 


2.1.3.7 


Control 




2,600 


650 


3,251 


2.1.3.8 


Tunnel 




512 


2,048 


2,560 


2.1.3 




Total 


29,116 


9,199 


38,316 
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Table B.6: Levels 4 and 5 of WBS and cost (in 2002 k$) of ME-section of SCL. 





SCL 




[\ /T O T~ t \ T 1 O 1 

IVlalcl lal 


T , ahnr 
JjaUUl 


Tntal 

lUtal 


2 14 1 


( ,vnn- A/f nniil p q 

\_y 1 U yJ 1V1 vJ\Jb LhvC-u 




3 837 


959 


4 796 






\z\ y\ Ire 


1 27Q 


320 


1 ^99 




2.1.3.1.2 


RF Cavities 


2 558 


639 


3 197 


2.1.4.2 


Refrigeration 




2,414 


603 


3,017 


2.1.4.3 


Warm Insertions 




925 


231 


1,156 




2.1.3.2.1 


Quadrupoles 


555 


139 


694 




2.1.3.2.2 


Steering 


92 


23 


116 




2.1.3.2.3 


BPM 


92 


23 


116 




2.1.3.2.4 


Vacuum Sys. 


185 


46 


231 


2.1.4.4 


RF Power Supply 




12,960 


3,240 


16,200 




2.1.3.3.1 


Klystrons 


3,888 


972 


4,860 




2.1.3.3.2 


Transmitters 


3,110 


778 


3,888 




2.1.3.3.3 


Circulators 


2,462 


616 


3,078 




2.1.3.3.4 


Controls 


1,814 


454 


2,268 




2.1.3.3.5 


Waveguides 


648 


162 


810 




2.1.3.3.6 


Couplers 


1,037 


259 


1,296 


2.1.4.5 


Electrical Distribut. 




350 


88 


438 


2.1.4.6 


Control 




2,049 


512 


2,561 


2.1.4.7 


Tunnel 




513 


2,051 


2,564 


2.1.4 




Total 


23,047 


7,685 


30,732 
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Table B.7: Levels 4 and 5 of WBS and cost (in 2002 k$) of HE-section of SCL. 





High-Energy SCL 




Material 


Labor 


Total 


2.1.5.1 


Cryo-Modules 




3,594 


898 


4,492 




2.1.3.1.1 


Tanks 


1,198 


299 


1,497 




2.1.3.1.2 


RF Cavities 


2,396 


599 


2,995 


2.1.5.2 


Refrigeration 




2,367 


592 


2,959 


2.1.5.3 


Warm Insertions 




832 


208 


1,040 




2.1.3.2.1 


Quadrupoles 


499 


125 


624 




2.1.3.2.2 


Steering 


83 


21 


104 




2.1.3.2.3 


BPM 


83 


21 


104 




2.1.3.2.4 


Vacuum Sys. 


166 


42 


208 


2.1.5.4 


RF Power Supply 




11,520 


2,880 


14,400 




2.1.3.3.1 


Klystrons 


3,456 


864 


4,320 




2.1.3.3.2 


Transmitters 


2,765 


691 


3,456 




2.1.3.3.3 


Circulators 


2,189 


547 


2,736 




2.1.3.3.4 


Controls 


1,613 


403 


2,016 




2.1.3.3.5 


Waveguides 


576 


144 


720 




2.1.3.3.6 


Couplers 


922 


230 


1,152 


2.1.5.5 


Electrical Distribut. 




330 


83 


413 


2.1.5.6 


Control 




1,864 


466 


2,330 


2.1.5.7 


Tunnel 




483 


1,934 


2,417 


2.1.5 




Total 


20,991 


7,060 


28,051 
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B.3 Conventional Facilities and Target Hill Cost 



Table B.8: Conventional facilities and target hill cost. 



Conventional Construction 


Quantity 


$ / unit 


Cost (k$) 


Upgrade U-Line Shielding 


275' 


5500 


1512 


Removals 


LS 




250 


Fill 


726350 cu M 


18 


13075 


Berm Stabilization 


LS 




800 


Primary Beam Tunnel & Liner 


1100' 


1200 


1320 


Decay Tunnel & Liner 


615' 


600 


369 


Target Area Footings 


600 cu yds 


500 


300 


Beam Dump Footings 


300 cu yds 


500 


150 


Retaining Walls 


800 cu yds 


600 


480 


Road Improvements 


LS 




500 


Target Area Shielding 


6200 tons 


350 


2170 


Power/high & low Voltage Distribution 


LS 




4250 


Water/Sanitary/Storm 


LS 




1350 


Power Supply Houses 








Lower Beam 40 x 100 


4000 sq ft 


250 


1000 


Upper Beam 20 x 80 


1600 sq ft 


250 


400 


Target Building 50' x 55' 


2750 sq ft 


500 


1375 


Near Detector Underground Facility 


LS 




6750 


Near Detector Building 30' x40' 


1200 


300 


360 


Cooling Tower 3.5 MW 


LS 




1850 


Cooling Tower .5 MW 


LS 




150 


Monitoring Wells 


5 


25000 


125 


Conventional Total 






38,536 
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Table B.9: Beam components costs. 



Beam Components 


Quantity 


q> / unit 


Cost (k$) 


39 Degrees Vertical Bend 


20 Dipoles 


100,000 


2000 


68.2 Degrees Horizontal Bend 


28 Dipoles 


1 A A AAA 

100,000 


2800 


Quads 


1 A 

14 


A A AAA 

90,000 


1260 


H&V Trim Magnets 


18 


45,000 


810 


BPM's 


12 


25,000 


300 


Current Transformers 


2 


35,000 


70 


Vacuum Components 


1200 


2000 


2400 


Trim PS's 


18 


25,000 


450 


Quad PS s 


14 


75,000 


1050 


Horiz. Bend PS 


1 


300,000 


300 


Vert. Bend PS's 


3 


150,000 


450 


PS & Instr. Controls 


300 Devices 


8,000 


2400 


Instrumentation -Flags/Loss Moni- 


25 Devices 


30,000 


750 


tors/etc 








Security Hardware 


3 Gates 


50,000 


150 


Horn Power Supply 


1 




3200 


Horns + spare hornl & Target 


LS 




538 


Components Total 






18,928 



Table B.10: Installation cost. 



Installation 


Quantity 


$ / unit 


Cost (k$) 


Installation Materials 


100 Devices 


5000 


500 


Remove/ Inst all Beam Dump 


25,000 tons 


20 


500 


Install Target Area Shielding 


6200 tons 


10 


62 


Magnet Installation 


77 


10,000 


770 


Vacuum Installation 


1200 ft 


200 


240 


Instr. Installation 


40 Devices 


10,000 


400 


Power & Tray Installation 




8 my 


800 


Security System Installation 




6 my 


600 


Horn Installation 




4 my 


400 


Installation Total 






4,272 
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Table B.ll: Conventional facilities direct costs. 



Conventional Construction 


Cost (k$) 


Conventional 


38,536 


Beam Components 


18,928 


Installation 


4,272 


Total direct costs 


61,736 
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